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The Interaction Rates of Stopped Negative Muons 
in Iron and Copper} 
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ABSTRACT 


An apparatus suitable for measuring the mean lifetimes of negatively- 
charged muons brought to rest in targets of Z<30 is described. Measure- 
ments made on iron and copper yield mean lifetimes of (204-+14) x 10-® sec 
and (157+ 20) x 10~ sec respectively. 


§ 1. IyrRopuctTIoNn 

Ix the ten years which have elapsed since the crucial experiment of 
Conversi et al. (1947) the interaction rates of negative muons have been 
investigated by many workers (see Sard and Crouch 1954). Recent work 
on parity non-conservation has again aroused interest in the muon and it 
seems likely that the interaction probabilities of negative muons will be 
of interest in determining the nature of the four fermion interaction and 
the structure of the nucleus. 

The apparatus to be described has been developed with a view to 
measuring the lifetimes of cosmic-ray muons in elements in the range 
Z<30, the limit being set by the channel width of the timing apparatus. 
Using a simple selection scheme it has been possible to determine the 
mean lives of negative muons in iron and copper by measuring the life- 
time distribution of a mixture of positive and negative cosmic-ray muons. 
The curve obtained is the sum of two exponentials of widely differing 
decay constants and it is possible to separate the two components with 
fair accuracy. 

§ 2. APPARATUS 


2.1. General 


Figure | shows the general arrangement of the apparatus. Incoming 
cosmic-ray muons enter the target by way of the Geiger counter tray G, 
and the scintillator S,. The disintegration products produce signals in 
S, while G, is in anticoincidence to reduce the background. In the 
present experiments, the target was 16 in. wide and 3 in. deep and rested 
on three pieces of angle iron. The iron target was circular in cross section 


and the copper target square. 


+ Communicated by Professor J. G. Wilson. 
¢ Now at Atomic Energy Research Establishment, Harwell. 
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The various parts of the apparatus were supported on a light frame- 
work which was, as far as possible, kept sufficiently far from the counters 
to render its effect on the counting rate negligible. The apparatus was 
assembled under a concrete roof, effective thickness 3 in., and no attempt 
was made at shielding from local radiation. 


S; 


Inches 


General arrangement of target and detectors. 


A negative muon coming to rest in the target, will produce signals in 
G, and $8, but not in G,. Neutrons and y-rays emitted downwards from 
the capturing nucleus will not produce signals in G, but may in S,. 
Neutrons can be detected in liquid scintillators by the ionization produced 
by knock-on protons, and y-rays by that from ‘Compton’ electrons. 
Where the negative meson decays, it is unlikely that the electron will 
escape from the target, but 8, can respond to bremsstrahlung from the 
electron. It is for this reason that the selection scheme does not reject 
events produced by positive muons coming to rest in the target, 


i 


Interaction Rates of Stopped Negative Muons Tit 


2.2. The Counters 


The trays G, and G, were made up of 18 and 17 Geiger counters (20th 
Century Electronics G60) respectively, the counters being connected in 
pairs or triplets to quenching units similar to that described by Hodson 
(1953). The outputs from the quenching units of a particular tray were 
connected to a mixer circuit which produced a standard pulse when any 
one, or more, of the counters fired. The counters were protected from 
daylight by thin aluminium trays and supplied from a source of positive 
K.H.T. Under these conditions, they were found to be quite reliable and 
the counting rates of the individual trays remained constant within the 
expected limits. 

The scintillation counter 8, was originally a liquid scintillator 
constructed in the following way. Two thin electrolytically polished 
aluminium sheets were cemented to a polished aluminium ring with Araldite 
D, cured for three hours at 60°c (see fig. 2). A glass window which allowed 


Scintillator 


a 


Schematic diagram of liquid scintillator Sj. 


the photomultiplier to ‘ see ’ the edges of the scintillator was cemented in 
at the same time. This container was filled with a solution of terphenyl 
and PopoP in toluene. It was cheap to make and very efficient and reliable 
in use. In about a year’s service, no evidence of ‘ poisoning ° of the scin- 
tillator was observed. When large plastic scintillators became commer- 
cially available, the liquid scintillator was replaced by a disc of Nuclear 
Enterprises (G.B.) Ltd. plastic scintillator Nel01. The disc, 16 in. 
diameter by 1 in. deep, was machined to a satin finish and recessed and 
polished in the centre to take a 2 in. photomultiplier. It was housed in 
a spun aluminium box and packed around with magnesium oxide. 

It was important that the detector 8, should be as thin as possible to 
reduce the numbers of muons stopping in it after being detected, since 


the lower portion of §, is effectively part of the target. Using such thin 


scintillators and a thin aluminium casing was effective in reducing such 
events to negligible proportions. The two scintillators described were 
tested for their detection efficiencies for muons. Using E.M.[. 6262A 


I2 
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photomultipliers, each detector was found to be more than 90% efficient 
for the detection of muons passing through its extreme edge, when the 
discriminator bias was set to produce a counting rate of 100 counts per 
second. 

The scintillator S, was désigned with a view to capturing neutrons 
(see Reines et al. 1954) and it is likely that a smaller disc-shaped scintillator 
would have been more efficient for the experiment described here. It was 
constructed from stainless steel, and was 12 in. in diameter and 12 in. 
deep. Before filling it was sand-blasted inside, and painted with a mixture 
of titanium oxide and water-glass. An electrolytically polished lid of 
thin aluminium was secured by a retaining ring and a lead washer. The 
scintillating liquid was viewed by three 6262A photomultipliers through 
glass windows cemented in with Araldite D. After a year’s use, the scin- 
tillator was inspected and found to be in excellent condition. The liquid 
was clean and free from any deposit. During the year, the coincidence 
counting rate remained sensibly constant. Some slight yellowing of the 
paint appeared to have occurred, but this evidently did not seriously affect 
the performance of the scintillator. It should be noted, that while the 
Araldite joints were resistant to toluene, the, were found to soften very 
rapidly when in contact with benzene or methyl! alcohol. 

The scintillating solution was prepared by bubbling nitrogen through 
toluene for severa] hours and then dissolving the terphenyl and PoPoP to 
concentrations of 3 g per litre and 0-1 g per litre respectively. Small scale 
tests showed that displacement of oxygen from the scintillator improved 
its performance by as much as 40%, in agreement with the results of 
Pringle et al. (1953). 

With 6262A photomultipliers, amplification of the output pulses was 
found to be unnecessary and, in fact, a crystal diode limiting circuit was 
incorporated in the photomultipler housings to limit the dynode pulses to 
10 volts. These dynode pulses were fed on to the grid of a CV2127 valve 
acting as a phase inverter and giving pulses of opposite polarity from anode 
and cathode. The positive pulse was fed directly to one grid of a pair 
of mixer valves and the negative pulse was fed to the other grid through a 
delay cable. The output of the mixer circuit was thus a negative pulse of 
duration equal to the delay introduced by the delay cable. This method 
of clipping the pulses was found to yield a square pulse free from any after 
effects. ; 

2.3. Electronic circuits 
General 

A block diagram of the apparatus is shown in fig. 3. Pulses from the 
anode of the S$; photomultiplier are fed through an inverting circuit to a 
fast trigger circuit, using EF P60 valves, similar to that described by Lewis 
and Wells (1954). The clipped pulses from the dynodes of 8, are fed to a 
diode coincidence circuit with a resolving time of 70 musec. If a coinci- 
dence is obtained, a gate is opened which allows a suitably delayed, 
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common pulse from the three anodes of the 8, photomultipliers to pass 
through to a second trigger circuit. The trigger circuits generate 
negative pulses 15 volts high and rising to maximum height in less than 
15 musec after the start of the initiating pulse. These pulses generate 
independently the front and back edges of the counter chronometer gate. 
As long as the counter chronometer gate is open, pulses generated by a 
20 Mc/s erystal-controlled oscillator are fed into a fast binary scaler. The 
counter chronometer is to be described in a separate paper. 


Fig. 3 


Delayed 
coincidence 
selector 


Final 
selector 


Mirrors 


okra 


Camera 
control 


Simplified diagram of electronic circuits. 
Q=Quenching circuit. C.F.=—Cathode follower. M.R.=Message register. 


The delayed coincidence selector and final selector circuit act as event 
selector circuits so that the camera circuit is only actuated when a G,S8,G, 
event is accompanied by a delayed S, pulse occurring within the gate 
length set by the delayed coincidence selector. The Geiger trays were 
necessary on account of the y-ray background from the walls of the 
building which produced a high rate of S, delayed S, events. Prompt 
coincidences were avoided by delaying the start of the delayed coin- 
cidence selector gate. The display of neon lamps which indicated the 
scalar count was viewed continuously by the camera (which had no 
shutter) by reflection in a mirror. When the camera circuit was fired, the 
mirror was moved slightly by a solenoid so that the image of the lamps 
was displaced on the film. The neon lamps were reset after two seconds 
and the camera wound on. In this way, any spurious counts of the scaler, 
or failures in the resetting circuits could be recognized when the film was 
read, by images produced on the wrong portion of the film before the 
mirror had moved. 
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The arrangement described above allowed the fast timing circuits to 
operate independently of the slower selection circuits. Thus every §,. 
pulse opened the counter chronometer gate and operated a paralysis 
circuit. The gate remained open for the duration of the trigger pulse 
unless previously closed by a pulse from S,. In either case, the infor- 
mation remained stored until the selection circuits decided whether to 
accept the event. If an event were rejected, the particular selection 
circuit which rejected the event triggered the resetting circuit which in 
turn sent resetting pulses to the counter chronometer scalers and also 
turned off the input paralysis. 


2.4, Haperimental Procedure 


When lifetime measurements were in progress, it was found that the 
rate of recording events was 25 per hour and it was necessary to change 
the camera film every 16 hours. Every eight hours, the output voltages 
of the power supplies were checked and the individual counter rates and 
the G,S, coincidence rate measured. A complete check of the counter 
chronometer using a pair of ganged pulse generators was made at the 
start and finish of every film. On one or two occasions it was found that 
the scaler had become unbalanced during arun. It was, however, always 
possible to identify the first false record and all subsequent records were 
neglected. 

Operation for about one week was necessary to obtain a lifetime 
measurement accurate to 10°% and it was easily possible to read a film 
carrying 400 events in one hour. 


2.5. Analysis 

During the run, the gate length of the delayed coincidence selector was 
set to 2-5usec. With this length of gate, the counts in the later channels 
are made up of positive muon events and chance coincidences, since the 
short lifetime of the negative muon ensures that they do not contribute 
significantly to counts in the later channels. The chance background 
was determined from the G,S8,G, rate and the known counting rate of 8,, 
and the positive muon contamination found by subtraction. It was then 
possible to extrapolate back to the early channels and finally obtain the 
numbers of counts in each channel caused by negative muons. The 
uncertainties introduced by this procedure in the final result were found 
to be small compared with the errors due to limited statistics. 

The results were calculated using the method suggested by Peierls (1935). 
In determining the root mean square error of the result, the signal to 
noise ratio (A in Peierls’ paper) was derived from the start of the first 
channel. As the positive muon contamination decreases exponentially, 
this is HO, to be an over-estimate yielding a slightly high value for the 
error. The signal-to-noise ratio of our apparatus found as above was 
about five. 

Meyer (1954) has considered some sources of systematic errors arising 
in a similar experimental arrangement and has shown that the soft 
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component of the cosmic radiation yielded spurious counts in the early 
channels of his detector. In view of this, we used no lead shielding, and 
also rejected the counts which occurred in the first three counter chrono- 
meter channels (i.e. 150 mysec), as inspection of Meyer's results showed 
that this eliminated the suspect counts. As the lower parts of S, and the 
upper parts of G, are effectively part of our target, it is necessary to 
consider the effect of their presence on the results we have obtained. 
Negative muons which stop in the light elements (C, N, H ete.) of the 
scintillators and counters will have a longer lifetime, in fact nearly 
2-2 usec, than those stopping in the target and will therefore appear to 
contribute to the positive muon background. 

The mass of the targets used was approximately 80 kg and the effective 
mass of aluminium and glass forming part of the target is estimated to be 
less than 1 kg. For the copper target, the effect of the iron framework of 
the apparatus and the target supports must be considered, an approximate 
calculation shows this to affect the result by less than 1%. We believe 
that systematic errors introduced by variations in the counter chrono- 
meter frequency and channel width are negligible. Observer errors in 
film reading are likely to be less with the display used than with a display 
on which measurement is necessary. Some films were read twice but no 
serious errors were detected. 

A check on the positive muon contribution was made in the first iron 
run when the gate was lengthened to 7-2 usec and the last 6-4 usec of the 
gate analysed. After subtraction of background, an exponential curve 
was obtained with a mean lifetime of (2-14--0-19) psec. By a rearrange- 
ment of the selection system, the mean lifetime of the natural mixture of 
positive and negative muons coming to rest in S, was also measured 
yielding the resuit (2-15--0-15) psec. 


§ 3. RESULTS 

The mean lifetimes of negative muons in iron and copper were found to 
be (204+14) musec and (157-20) musec respectively. Figures 4 and 5 
show the differential decay curves obtained. The table shows some of 
the results obtained by workers in this field and also the interaction rates 
assuming that the decay probability of negative muons is (2-22 psec)~', 
Bell and Hincks (1951). 

It is clearly seen in the table that there are discrepancies in the results 
from different laboratories, which cannot easily be explained in terms of 
the quoted errors. A few remarks on this subject may be useful here. 
Meyer (1954) has objected to some of the results of Keuffel e¢ al. (1952) on 
the grounds that their apparatus was troubled by spurious counts in the 
early channels of their timing apparatus. This led to a decrease in the 
measured lifetime and thus to high capture rates. He therefore repeated 
some of their measurements and for this reason, his result for copper is 
presumably more reliable than theirs. Meyer did not repeat their 
measurement on iron, but a similar fault in their apparatus would explain 
the discrepancy between their result and the last three results in the table. 
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The Mean Lifetimes and Interaction Rates of Negative Muons (Captured 
in Iron and Copper 


Mean lifetime | Interaction rate 


Element (x 10° sec) (x 10-* seo-) . Reference 
Fe 26 210+60 4-31+1-4 Benade (1953). 
163-27 5-68+1-0 Keuffel e¢ al. (1952). 
242415 3-68-40-23 Barrett et al. (1957) 
192-20 4:76-0:54 Lederman and Weinrich (1956) 
{202+7 4-49-+0-15 Sens e¢ al. (1957) 
§204-+ 14 4-45-+40-33 Present work 
Cu 29 116+9 8-17-L0-67 Keuffel e¢ al. (1952) 
172+8 5:36+0-27 Meyer (1954) 
7157412 5:9240-49 Lederman and Weinrich (1956) 
t138+5 6:79-0:22 Sens et al. (1957) 
§157-+20 5-92-+0-80 Present work 


+ Assuming the mean lifetime of negative muons in carbon to be 
(1-92 -+.0-04) psec. 
* Assuming the decay probability of negative muons to be 
[(2-21 + 0-02) x 10-®]-? sec“!. 
§ See note added in proof. 
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The recent results of Sens ¢/ a/. (1957) are most interesting as repre- 
senting a series of measurements on a large number of elements. It is 
noticeable that their measurements on copper yield a decidedly higher 
result for the interaction probability than other workers (if the Keuffel 
result is ignored for the reasons stated above). The shortness of the 
publication by Sens e¢ a/. makes it difficult to comment on this discrepancy, 
but it is interesting that in their comparison of their experimental results 
with a theory by Primakoff the fit is very good throughout the range of 
measured elements except in the case of copper where the predicted capture 
rate is 5-45 x 10® sec™?. 

The value of the interaction rate in iron (4:45-++0-26) x 108 see! taken 
in conjunction with that of Weinrich and Lederman for the interaction 
rate in calcium yields a value for the ratio of the interaction rates in iron 
and calcium of 1-77--0-177. A recent paper by Tolhoek and Luyten 
(1957) suggests that this ratio depends on the nature of the interaction 
between the four fermions involved in muon capture. Tolhoek and 
Luyten predict ratios of 1-08 and 1-70 for the scalar—vector interaction 
and tensor—axial vector interaction respectively using a potential well of 
radius 1-410 At?em. It would therefore seem that experimental 
results weigh against the pure scalar—vector interaction—a conclusion 
which is supported by the results of Sens et al. in other elements. Further 
measurements of other elements are in progress. 


Note added in proof.—By making use of more channels when calculating 
the positive muon contribution to the counts, we have been able to reduce 
the uncertainties in the estimations of the positive muon backgrounds and 
thus in the final results. As the extra channels used contained counts 
caused by negative muons, the estimation of the numbers of positive muon 
events was made by successive approximations. Calculated in this way, 
the mean lifetimes of negative muons in iron and copper were found to be 
(201+11) mysec and (146--17) musec respectively. 
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of a Classical Plasma} 


By 8. F. Epwarps 
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ABSTRACT 


A variational calculation is given for the free energy of a fully ionized 
classical plasma. The method avoids redundant variables and allows for 
cross terms between collective and individual modes of oscillation, though 
these turn out to be small. A relation is given for the shielding distance in 
terms of the Debye—Hickel radius and the interparticle distance, it being 
of the order of the former. 


§ 1. INTRODUCTION 


THE partition function of a classical plasma cannot be calculated in a 
straightforward way because the infinite range of the coulomb potential 
will not allow any expansion of the exponent in the formula 
Z=) exp (—AH/kT). 

This problem is understood physically in that the gas shows two kinds 
of behaviour: individual particle motion, and collective, plasma oscilla- 
tions. It is not easy to consider the collective oscillations as built 
up from individual particle motion. In many ways the two motions 
appear almost independent of one another, but of course the gas has 
only 6N degrees of freedom (N electrons, NV protons) and one cannot 
introduce more coordinates than this. However the concept of collective 
coordinates is only a means to the end of calculating physical quantities, 
and there is no reason why the same division of the degrees of freedom 
between collective and individual particle coordinates should be used 
throughout a calculation. In this paper the coulomb potential will 
be separated into long and short range potentials using a parameter 
h which is a variable of the calculation. The short range part can then 
be treated by expansion into clusters, a particular division between 
collective and individual coordinates being used for each cluster. The 
final answer should not depend upon the exact manner in which the 
separation takes place and so is stationary against a small change in h, 
which is used to find h. It is found, as is to be expected, of the order of 
the Debye—Hiickel radius. This calculation is much simpler than that 
of plasma excitation, since the velocities are independent variables in 


+ Communicated by Professor R. E. Peierls, F.R.S. Work partly performed 
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the partition function, and reliable estimates of the various approxima- 
tions customary in plasma theory can be made. They are all good. 
The contributions to the free energy from collective and individual 
motion are comparable. In order that this calculation can be pursued 
to higher approximation, quantum effects must be invoked to avoid 
the electrons and protons coalescing. These effects are not treated 
in this paper as they are negligible in the fully ionized gas discharge, 
but could be included if required. 


§ 2. Tur EVALUATION OF THE FREE ENERGY 


The free energy of the system is defined as 


—F=kxT log Z a le Se Ae 
where Z is the partition function 
Z| uS | Hare, d3z,exp —H(l1, 2,..., N)/KT’. (2) 
For a system of N protons and N electrons 
H=)> 4mv, Pay: 4Mv 2+ 2s Ticetak te) ee 
x%6 


?, J referring to electrons, x, 8 to protons. The velocities integrate 
out straight away and contribute additively to #. Henceforward F 
will be taken to refer to the remaining potential energy part. 

The simplest separation into long and short range forces is to write 


EoeeOxD (=F) 1 exp (—7/h) 
: ae Strela tn Sear (4) 
If the latter is fourier a in a box of side R, it gives 
1 exp (—7/h) 
- a aa = PS exp ( (2m7k . r/k){k —(k?-+(R/ [2mh)?}-4 (5) 
if 
N 
Pe= > exp (27k . r;/R) 
‘x (6) 
Ti SAGED (2atk . r,/R 
and 
(R/2h)2k-*(k? + (B/2h)?)-1= =) iets, cee) 
then 
1 wk r,;[R) a 
x{e- Sed See >>I ed 


It is convenient to use cate N of the p's (or subsequently N-6 etc.), 
so define a potential 4(r) 


1 
> H(ri5) is) aye Dd(ruPr™ —N)¥(k) ‘ . : . . ° (9) 


ae) 
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where the sum now runs over N vectors k, i.e. inside a sphere of radius 
(3N/4)'3 in k space. The part corresponding to the sum outside this 
sphere to infinity corresponds roughly to a potential which is significant 
only inside the interparticle distance, and should be considered as part 
of the short range potential. Thus now 
1 
ee) + (5-40 ) =¢$(r)+x(7) eae (10) 
is used as separation. 
The potential energy now becomes 


»> (xis Xop— 2X ia) + DP(k)(pupn*—N +o,0,.* —N 


a#B 
— Py Fy — p,*o,)- eth cee eee StL) 
If p,—o,=p,, the latter becomes 
YVP(k)( Per —2N). ee ee LS) 


Now consider the calculation of the free energy when firstly the latter 
and then the former terms are neglected. The former can be treated by 
expansion according to the formula 

<exp (—H,/xT))=)>(—IeT)™(H,” > pean ee SATs) 
which can be rearranged into 

exp{—H,,/«T+3((H,")— (H, )?)/e°T?— 3((H,®) —3(H, )<H,*) 
+2¢H,)*)/KP1® (14) 

so that 

F=(H,)—3}((H2)—(H, eT? 6... (18) 
and for H,=e?>(xi;+Xag—2Xia) this series converges well provided that 
h>ro, A>ro where r,°>N=R3, N=xTr,3/4re7. It gives 


(a) 0) 
(A 2 )=Nr,Ph/(87A*) ee here fice ECL} 
Gla Bie east Ia oe Os & Sieur ce we Ala en ee aM Chih) 


where ¢ involves the quantal cut-off required to avoid an ultra-violet 
catastrophe. However assuming the / will be of the order of A, and 
physically the value of r,A~1<1 this term is already very small, and the 
expansion will be pursued no further. The change in separation of 
long and short range potentials from exp (—r/h)r-' to d(r) makes a 
change of order (r,/h)? in the free energy and confirms that exp (—r/h)/r 
can to all intents and purposes be taken as the shielded potential. 
If the collective part of the Hamiltonian is now considered, it gives 


[...] Mae, exp (—L(rare* 20) ¥ (Hk). pers) 


If the p,’s are now taken as the basic set of coordinates (with the 
other combination p,+o,, which does not enter the Hamiltonian) a 
Jacobian is required from the 2’s to the p’s. Although this is well 
known to be exp—{p*p/2N} (e.g. Bohm and Pines 1952), a derivation 
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is given in an Appendix which also estimates the errors involved. 
Strictly the |p, | cannot exceed 2N, but the integral (18) is strongly 


convergent and can be taken to be 
ee IT dp, exp{ —Dy byt (F/T + (2N)-)+2NY eT } 
kin sphere 
=exp [NPT r(P/ eT + (2N)- 2) (KT XP. og on, Oo 


The contribution to the free energy is then, omitting additive constants 
—F 

Now consider the whole problem. The series expansion cannot now 
be simply applied with the above evaluation of the collective contri- 


bution except in the first term, which is just the collective contribution 
alone. In the term (H,?) integrals are in the form 


[oriy) Px, de, 
collective term, and one is using all the coordinates already in the collective 


term. This is easily overcome by separating out the 7,9 coordinates, 
say electrons 1 and 2 without loss of generality and redefining 


long 


ages exp (277k . SAE bs exp (277k . r,/R) 3 (aa 

Then : 

Dd (PP — 2N yb (k) > 

Dd {oi Dy *— py, (exp (— 27k . r,/R)+-exp (—27ik . r,/R))+c.c. 
t+exp {277k . (r;—r,)/R}+e.0.—2N+2} . . (22) 
The p,’ can now be integrated away, and the only final change lies in the 
2N—2 rather than 2N is the exponential, and that the & sphere now 
should have a radius of order (N—1)'3 rather than N'/3. Both of these 


are order N~! down from the principal terms, and in effect the two inde- 
pendent calculations can just be added together: 


—F=N>P+kKTY log (W/ceT+(2N)-1)+(Nry2h)/(87A4eT) 2. (238) 


§ 3. THE DETERMINATION OF THE SHIELDING RaprIus 


The free energy mainly comes from short range encounters of the 
particles and from long wavelength fluctuations of the charge density. 
Between them lies a region which roughly speaking can be treated in either 
way and somewhere in this region the diversion is made which involves 
the radius h. Ifh is increased by 65h, this amounts to removing some part 
of the short range potential and including its effect in the long range 
potential. It is reasonable that there is a range of values of h for which 
it makes no difference whether this part of the potential is treated in one 
way or the other. It has been checked that the actual F(h) has a broad 
minimum and the value of h giving this minimum represents the best 
value of h, i.e. that which makes the calculation most accurate. It has 
been assumed that taking more terms into the calculation will not 


a 


So 
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materially alter this position, and this certainly appears to be the case. 
As it stands, however, the method does not have the power of a minimal 
principle. 
After algebraic manipulation the equation d//dh=0 becomes 
19° *h-1= 167 > (k?[ 2h R]P+ 1) -2{N 4-422 Rr 3h? k2( Qh) R)2+ 1-3}. 
(24) 
The sum over k can be replaced by the integral inside the sphere 
(4/3)7k>max=N, and putting k=k’(4/3)1/3(27)-1, (24) becomes 
SrAZzA 2a(32/4)1/8 ; 
1= ries (kPh? rg? 1)-1(1-++ k?A279-2[k2h27 2+ 1])- 1k? dk. (25) 
Provided that A >r9, which is the case under consideration, the upper 
limit of this integral can be extended to infinity and finally the result is 
LANIER WY pasta ie salen 
(S)=S[ protect ea . 26) 
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which gives 2h?~A?, i.e. the shielding distance is of the order of the 
Debye—Hiickel radius, and the free energy is now given by (23) using h 
from (26). 

It is interesting to note that the contribution to the free energy from 
both types of contribution is of the same order: (V/A°)kT' (V the total 
volume). 


§ 4. CONCLUSION 


It appears from this calculation that the usual approximations of 
plasma theory are justified in thermal equilibrium, and that the method 
given here gives a reliable calculation of the thermodynamic quantities of 
a fully ionized gas, in the region where quantal effects play no part. The 
problem of the correlation energy can be approached in a very similar 
manner except that since the free energy is a number, so also is h, whereas 
the correlation energy, being a function of r requires an A(r). It is hoped 
to publish a treatment of this in due course. The other problem remaining 
in thermal equilibrium is that of the correlations in space and time of the 
charge density in an equilibrium state, and the frequency of oscillations 
associated. It is hoped to treat this also. 


SPN DX 


It is required to find the Jacobian of the transform «—p. For 
simplicity consider the one dimensional transformation 


N/2 
Py== > exp{2atk . x,/R). 
N/2 
The Jacobian is the ratio of the elements of volume dz,...... dx to 


dp_njg-++-> dpyj. Now this can be regarded as the probability of 
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finding p_yj-.+++ Py When the 2,..... ay is random since if P is the 
probability of a distribution 
Pla, ...@y) de,...dty=dp_yo...Apyel(p_ wie --- Pm) 


and if the 2’s are random P(«,...«y) is unity. This probability is the 


~ P_N/2 PA/2 : 
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The value of this integral lies almost all near v0 and the x integration 
can be performed by expansion and recombination in the exponential. 
This is essentially the method of stationary phase. Then the first integral 


1s 
p00 +00 N : 
Hem {" ny] exp| — Bie — = Db 
Ni meee 
Wh Wohoda+d-+e+d) | =IT exp (—p,p,*/N) @ 
th 
when JN is large this becomes 
(7/N)* exp (— >'p,p;,*/N) eT rs hi 


Note that the form (27) will vanish when the ’s lie outside the interval 
(— 3K, 3) and so P(p) vanishes over a certain part of the p space. This 
is represented by the terms in the expansion, which are so small that 
effectively the integral can go to infinity. 
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ABSTRACT 


Tn the face-centred cubic lattice a stacking fault, caused by the condensation 
of vacancies on a {111} plane, is converted into a dislocation source of 
Frank—Read type if it is swept by a half-dislocation of a kind that restores 
lattice order, Dislocation sources will also arise in other major lattice types 
through the condensation of vacancies on close-packed lattice planes. 
It is proposed that these are the dislocation sources which release a major 
proportion, if not all, of the dislocations formed during plastic deformation. 
Besides several minor advantages, the introduction of this kind of dislocation 
source can account for all the important glide systems found in erystals of 
the most common lattice types. 


§ 1. [InTRODUCTION 


For a number of years now the theory that the dislocations formed 
during plastic deformation are released from a random network of 
dislocations (Mott 1952) in which the links act as Frank—Read sources 
of dislocations (Frank and Read 1950) has been accepted almost 
universally. This is not surprising since this theory is at the same time 
simple and theoretically convincing. Not only are there theoretical 
reasons to account for dislocation networks, even in the most carefully 
prepared specimens (Frank 1950), but moreover Hedges and Mitchell 
(1953) found evidence of dislocation networks in silver bromide crystals. 
It seems quite certain that, given suitable stresses, such a network can 
release dislocations by the Frank—Read mechanism. An important 
question, however, still remains to be answered, namely whether in 
actual fact all, or at least most of, the dislocations produced during 
plastic deformation are formed in this way. 

Since the above theory of three-dimensional dislocation networks was 
proposed the number of published experimental observations which appear 
to be in more or less severe disagreement with it has increased. In the 
following five such observations are listed. 

(i) At small strains, before slip bands appear, face-centred cubic pure 
metals show surface striations, usually not fully resolved even in the 
best electron micrographs. These have been named ‘fine slip’ 
a a aT 
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(Brown and Honeycombe 1951) or ‘elementary structure’ (Wilsdorf 
and Kuhlmann-Wilsdorf 1951). Measurements of this structure on 
super-purity aluminium crystals suggest—for every one of the twelve 
possible glide systems—the presence of about 107 sources/cm?, which are 
capable of acting at stresses below 1 kg/mm? (Wilsdorf and Kuhlmann- 
Wilsdorf 1952, Kuhlmann-Wilsdorf and Wilsdorf 1953). Such a source 
density is rather too high to be easily credible. 

(ii) Recently Toaka and Sakata (1957) observed that strained specimens 
of Cu,Au alloy exhibit the elementary structure when deformed while 
in the ordered state, but show the long sharp lines first described for 
a-brass (KuhImann-Wilsdorf and Wilsdorf 1953) when disordered. The 
latter type of surface structure would indicate a very much smaller 
number of active Frank—Read sources than the former, and it is difficult 
to see how a simple heat treatment could have such a profound effect 
on the structure of a dislocation network. 

(iii) A similar unexpected effect of heat treatment on surface structure 
is observed on pure aluminium and copper. These metals show elementary 
structure and slip bands when well annealed (Kuhlmann-Wilsdorf and — 
Wilsdorf 1953) but produce a slip line structure resembling that of 
x-brass when deformed after rapid quenching from a high temperature 
and subsequent ageing (Maddin and Cottrell 1955, Maddin 1957 (private 
communication) ). 

(iv) In a random network of dislocations a considerable fraction of the 
links, acting as Frank—Read sources, would give rise to the so-called 
“homogeneous siip’ (Suzuki and Fujita 1954, Suzuki 1954, Bilby 1955), 
but a thorough study of the surface structure on «-brass (Wilsdorf and 
Fourie 1956) failed to reveal any homogeneous slip. 

(v) Gilman and Johnston (1956, 1957) Gilman (1957), in a series of ingen- 
ious and careful experiments, found no evidence of pre-existing dislocations 
at the centre of newly-formed dislocation loops, and concluded that, in 
crystals of Lif’, dislocations are not released from (network) Frank—Read 
sources. 


§ 2. Srackine Favuuts as DisLocaTion SOURCES 


A dislocation represents a very strong concentration of elastic energy. 
This fact makes it unlikely that dislocations can arise through the action 
of stresses comparable to the typical yield stresses of pure face-centred 
cubic metal crystals, and, because a dislocation in order to be stable 
at such stresses has to be about 10~4 cm long, thermal fluctuations cannot 
materially assist in their formation. Further, since the most carefully 
prepared specimens usually have the lowest yield points, local stress 
concentrations, such as might exist in the vicinity of sub-boundaries, 
impurities and other lattice defects, cannot be responsible for the 
formation of dislocations in metal crystals at low stresses. 

The obvious remaining way, then, to explain the profuse generation 
of dislocations taking place at low stresses in metals is to seek their origin 
in some kind of lattice defect which can be assumed to be present in 
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large numbers in every real crystal. This was done most successfully 
by Frank and Read in that mechanism which carries their name, after 
a previously-developed model of dislocation multiplication by reflection 
(Frank 1948) had for some time been widely accepted as the correct 
explanation for the origin of dislocations. . 

Dislocation reflection may well take place under certain conditions, but 
it is not now believed to be the normal process of dislocation production 
in ductile crystals, mainly because calculations showed the required 
stresses to be too high (see Cottrell 1953). Further, the predicted rapid 
formation of slip lines, and the corresponding sudden strain increments 
could not be verified experimentally (Leibfried 1950). Another possible 
model for the generation of dislocations, explored by Frank and Stroh 
(1952), is the formation of dislocations at the ends of kink bands; but 
this requires kink bands of a strength which cannot possibly be expected 
in pure, undeformed crystals. 

Since Seitz (1952) studied the interactions between vacancies and 
dislocations much evidence has become available to indicate the profound 
influence which vacancies—singly, in groups or condensed into stacking 
faults—can have on the plastic properties of metals (see the review article 
by Cottrell 1957). Therefore it seemed promising to search for a process 
by which vacancies would be responsible for the creation of dislocations. 
The simplest way for this to happen in the face-centred cubic system 
would be by the condensation of vacancies into two adjacent, continuous 
and similar layers in the {110} plane, bounded at one end by a (110) 
direction. This would yield a normal edge dislocation (Barnes 1954) or, 
more accurately, an edge dislocation forming part of a prismatic dis- 
location (Seitz 1952). However, there is no apparent reason for vacancies 
to condense in this peculiar fashion, but instead it is highly likely that 
they will condense in {111} planes, thereby forming stacking faults 
bordered by Frank sessile dislocations (Read 1953). 

Seitz (1952) believes such vacancy condensation to be possibly * a very 
important means by which dislocation rings are generated in practice, 
for example, during cooling of a melt in which there is a high density of 
vacancies...’. Again, this idea of dislocation formation through the 
condensation of vacancies has been employed to account for the origin 
of sub-grain boundaries (Teghtsoonian and Chalmers 1951, 1952, Frank 
1956). Moreover, the same kind of stacking faults could, and probably 
do, arise through growth accidents during crystallization from the melt 
or during recrystallization. Hence it is reasonable to assume that 
significant proportions of stacking faults are formed in every real crystal. 
although probably not enough to be detectable by x-rays (Smallman 
and Westmacott 1957). Further it appears to be most likely that the 
stacking faults in face-centred cubic crystals will tend to be bounded 
by (110) directions, these being the close-packed lattice directions. 

It is a property of the stacking faults under discussion that order is 
restored if the area of the fault is swept by a half-dislocation having any 
of three possible slip vectors in (211) directions. (‘The other three, 
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of the total of six slip directions—-counting forward and reverse motion 
separately—which lie in any {111} plane are excluded because they would 
cause the atoms in the two adjacent layers to face each other.) Any such 
half-dislocation will combine with the Frank sessile dislocation according 
to 4[211]+4{111J=4[011], thus forming a total dislocation of normal 
type (Read 1953). Depending on the slip direction of the half-dislocation, 
the resultant dislocation may have any of the three slip directions of 
{011) type which do not lie in the plane of the fault (figure). From now 
on a dislocation formed in the manner described, and encircling the site 
of an eliminated stacking fault, will be called an R-dislocation—for 
‘resultant ’ dislocation. 


The geometry of slip planes, and slip directions of complete and 
of half dislocations in the face-centred cubic lattice. 


As was argued before, stacking faults in face-centred cubic lattices will 
tend to be bounded by ¢110) directions, and hence the R-dislocations 
will mostly consist of straight segments parallel to the three (110) 
directions which lie in the plane of the original stacking faults. These 
directions, together with the slip direction of the R-dislocation, define 
the possible slip planes on which the R-dislocation could glide. From 
the figure it can be seen that only two of the three (110) directions in the 
plane of any particular fault share a {111} plane with the slip direction, 
i.e. share that plane which experimentally is known to act as slip plane. 
Consequently R-dislocations can act as Frank—Read sources, such that 
those parts which lie in a favourably oriented {111} plane can form 
expanding loops. Those parts which are parallel to the third, unfavourable 
{110) direction or which lie in a little-stressed {111} plane provide the 
immobile or dragging dislocation nodes, 
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§3. Some Numericau Estimates ON R-DIsLocaTIONsS 

in the preceding section we saw that R-dislocations will arise wherever 
half-dislocations are formed in a stacking fault, and it should now be 
investigated under what conditions they can be expected to form. 
Obviously the creation of half-dislocations in stacking faults is a much 
easier process than the formation of glide dislocations in an undisturbed 
lattice. This is so because a half-dislocation has a considerably smaller 
energy and because it is strongly attracted towards the still-faulted 
regions, so that it moves as if under the influence of a strong shear stress 
tending to increase the diameter of a newly-created loop. The resolved 
shear stress corresponding to this attractive force depends on y, the 
stacking fault energy. In face-centred cubic metals its magnitude is 
estimated between about 20 ergs/em? and 200 ergs/em?, whereby the 
lower figure would be applicable to copper, the higher to aluminium 
(Seeger and Schéck 1953, Smallman and Westmacott 1957). The 
equivalent shear stress of a stacking fault is easily calculated as follows. 

If a lem dislocation line with slip vector b moves under the 
action of a shear stress 7 through the distance s the mechanical 
work done thereby equals tsb. The energy released if the same 
dislocation moves over a stacking fault, leaving order in its path, is 
given by ys. Equating the two expressions, the equivalent shear 
stress becomes 7,—y/b. With b=2x10-$cm and y=20 ergs/cm? this 
yields 7,=10° dynes/em?=10 kg/mm?, and with y—200ergs/cm? the 
equivalent stress is 7,—100kg/mm?. The typical yield stress for pure 
face-centred cubic metals is about 0-1kg/mm?. According to a 
calculation and curves given by Frank (1950) it appears that a stacking 
fault energy of 200 ergs/cm? will lead to the spontaneous nucleation 
of half dislocations, but not a fault energy of 20 ergs/cm?. 

An additional factor which may have to be taken into consideration 
is a possible interaction between the sessile dislocations and the half 
dislocations. Since the resulting undissociated dislocation has about as 
much energy as the two single dislocations together before the reaction, 
there is no first-order effect. A second-order effect may exist, though, 
as the resulting dislocation—where it lies in a {111} plane—will dissociate 
into half-dislocations of type $(211), thereby reducing its energy ; and 
also because of changes in the core energy. Superficially one might think 
that this interaction could assist in the formation and initial growth of 
a half-dislocation if it were nucleated at such a small distance from the 
sessile dislocation that parts of them could combine while the half- 
dislocation was still very small. 

Considering all these points together it seems probable that in those 
cases in which half-dislocations are not formed spontaneously great 
numbers will be formed through thermal activation during the unavoidable 
periods of heat treatment. 

For specimens still containing a majority of unconverted stacking faults 
it is to be expected that half-dislocations will be formed through the 
action of the stress fields of dislocations which either pass through the 
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faults or move past them at close distances. From the first few dis- 
locations a chain reaction will result since the newly-formed R-dislocations 
can release loops which in turn will trigger more R-dislocations and so on. 

The process outlined is possible because the stresses near the centre 
of a dislocation reach values, sufficient for the formation of half 
dislocations in stacking faults, within distances of 10 atomic diameters 
or more from the dislocation axis, and because the critical radius of a 
half-dislocation in a fault is quite small. It remains to be seen what 
types of dislocations are likely to result in this way. 

As was found earlier the slip direction of an R-dislocation cannot le 
in the plane of the original fault. Therefore, if a dislocation approaches 
a stacking fault in its own slip plane or cross slip plane the resulting 
R-dislocation cannot have the slip vector of that triggering dislocation. 
Since normally an R-dislocation can release dislocations in both planes 
which contain its Burgers vector we see that, at the very minimum, 
sources on four different glide systems must be produced. Moreover, 
remembering that the active stresses strongly depend on the relative 
position of fault and triggering dislocation, on its type—whether edge, 
serew or mixed—and possibly on its speed (which may determine whether 
the nucleation of the half-dislocation takes place ahead of the approaching 
dislocation or behind it), one comes to the conclusion that the passage of 
just one dislocation loop through a random array of stacking faults can 
lead to R-dislocations on all possible slip systems. 

Another problem to be considered concerns the probable size and 
number of stacking faults in a crystal. An estimate of the total faulted 
area can be made from the concentration of vacancies at the melting 
point. This is approximately given by c > exp (—U/kT) with U the 
energy of formation of a vacancy. If all the vacancies which are present 
in a crystal at its melting pomt would condense into stacking faults the 
total fault area per cubic centimetre could be found as 


A = (1/d)[exp (—U/kT)] 


where d is the distance between closely-packed lattice planes. This 
formula will give a somewhat too high value because a fraction of the 
vacancies will not find their way to stacking faults but will be absorbed 
by other types of sinks for vacancies. On the other hand there will be 
an additional—and so far unknown—contribution due to faults arising 
from growth accidents. Thus the above formula seems to be adequate, 
particularly since there exists some uncertainty about the value of U 
for most crystals, which probably introduces a larger error than that 
due to vacancies not captured by stacking faults. 

To make an estimate as to the average size of the stacking faults is 
a difficult task. It is dependent on a number of parameters, as for 
instance, the thermal history of the specimen under consideration, the 
stacking fault energy, the activation energy for the nucleation of a stacking 
fault, the binding energy between vacancies to form groups, and the 
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diffusion constant of single vacancies and aggregates of vacancies. 
However, it is still possible to give limits and probable trends: the 
diameters of stacking faults probably range from several microns to 
less than 50 A. An upper limit’ of stacking fault diameter is to be expected 
because the faults will hardly become larger than the cross section of the 
sub-grains in their parent crystal. A lower limit is necessary because 
below a critical size the energy of all contributing vacancies together 
becomes smaller than the energy of the stacking fault plus the energy 
of the dislocation surrounding it. 

In order to estimate the lower limit of stacking fault diameter we 
compare the surface energy of a disc-shaped hole—the potential stacking 
fault before its sides have collapsed—with the energy of the fault plus 
surrounding dislocation. The energy of a dislocation loop of very small 
diameter, D, is about equal to Gb?D (Frank 1950), with G the modulus 
of rigidity, and the energy of the stacking fault is equal to 47D*®y. The 
sum of these two energies must be equal to, or smaller than, $7D®c 
where o is the surface energy of the metal considered. From this 
we get (D/b) = 4Gb/[7(20—y)]. For aluminium this formula gives 
D=104A. 

The wide variations in size and number of R-dislocations, which, 
according to the preceding argument, are obviously possible, can account 
for the observed great differences in the plastic properties of specimens of 
different metals, or of specimens of the same material but with a different 
history. For instance one will expect that in order to generate large 
stacking faults long periods of heat treatment will be necessary, and that 
specimens rapidly cooled from the melt will tend to have many small 
faults. So qualitatively one sees why single crystals should have low 
yield points (slow cooling leads to small supersaturations, fewer but 
larger stacking faults, yielding R-dislocations able to act at low stresses) 
while rapidly solidified specimens should be much harder. 

As was argued above, the longest segments of R-dislocations should be 
comparable in length to the diameter of the subgrains. Approximately, 
then, the yield stress 7, should be connected with the average diameter 
of the subgrains, L, as 7,=Gb/L, where G is the modulus of rigidity 
and 6 the Burgers vector. This formula was previously derived by 
Bragg (1942) in a quite different way and found to be in satisfactory 
agreement with experimental data. 


§ 4. Fine Sire, Suip Banps, «-BRass STRUCTURE AND CROSS SLIP 


Qualitatively the interpretation of the phenomena of elementary 
structure and slip bands will hardly be affected by the introduction of 
a new type of dislocation source. A theory developed earlier (Kuhlmann- 
Wilsdorf et al. 1952) seems to be adequate. There it was argued that at 
stresses near the yield point a multitude of dislocation sources start 
to act almost simultaneously, and that the dislocation loops emitted by 


182 Doris Kuhlmann-Wilsdorf on the 


them spread out fairly slowly until they are stopped through the inter- 
action with loop sequences from neighbouring sources. Because of the 
high density of active sources the mean free path of the dislocations 
remains so small that only comparatively few loops, say up to 30, can 
be emitted from any one source, leading to the observed shallow, 
closely-spaced elementary lines. 

In the same paper the formation of slip bands was explained as follows : 
the first few dislocation loops from any elementary lines which, acciden- 
tally, have become particularly strong can first encircle, and then force 
back into their respective sources those dislocation loops which block 
their path. The annihilation ‘of the blocking loops aids the action of 
other, suitably situated dislocation sources which in turn allows their 
loop sequences to grow strong enough to repeat the process, and so on. 

The interpretation of the «-brass structure, which we believe appears 
on alloys with a short range order, has already been given in a recent 
lecture (Kuhlmann-Wilsdorf et al. 1957). The significant difference 
between metals with an elementary structure and metals with an «-brass 
structure apparently does not lie in the number, quality or distribution 
of the dislocation sources, but in the changes caused on a slip plane by 
the passage of a dislocation. In an alloy containing a short range order 
- the yield stress is higher than the stress which would be necessary to 
activate the available dislocation sources in a disordered lattice. This 
is so because, by their motion, dislocations disturb the order along their 
slip planes, and thus the external stress has to provide this ordering 
energy in addition to that needed to move the dislocations through a 
disordered lattice (Fisher 1953, Cottrell 1953.a, Logie 1957). Oncee slip 
has started on any lattice plane the disordering proceeds at a decreasing 
rate, and complete disorder is reached after the slip on a plane is about 
equal in magnitude to the diameter of the average domains. Conse- 
quently the stress necessary to move a dislocation along a slip plane is 
highest for the first one to cross that plane and then diminishes until 
the stress corresponds to complete disorder. Therefore, although an 
alloy exhibiting the «-brass structure probably contains stacking faults 
of similar sizes and numbers as a pure metal—which may or may not 
have already been converted into R-dislocations—slip will only start 
at the threshold stress which is necessary to destroy the short range 
order ; but once slip has started from a source it will be propagated 
at a high speed over comparatively large areas of the slip plane, thus 
causing the long and sharply defined slip lines typical of «-brass. 

‘Clear’ cross slip, i.e. sharp slip lines abruptly changing from one 
slip plane to the corresponding cross slip plane, is prominent wherever 
the «-brass structure is observed, but, according to electron micrographs, 
it is exceedingly rare in pure metals (Wilsdorf and Kuhlmann-Wilsdorf 
1952, Noggle and Koehler 1957). This can be easily understood on the 
basis of our argument that sequences of dislocations in pure metals are 
obstructed by slip from neighbouring sources even at low stresses, 
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whereas in «-brass and similar alloys the sources which act at small 
strains are unimpeded by other sources. 

The straight segments of the R-dislocations which act as Frank—Read 
sources are, in face-centred cubic metals, mostly linked to one segment 
on the corresponding cross slip plane. Where there is little freedom of 
movement for the dislocations, as in pure metals, this will be of little 
consequence, but wherever the R-dislocations emit loops, unhindered by 
neighbouring dislocations, there will be a probability that the dislocation 
loops spreading out on the active slip system will, at the same time, 
drag out loops on the cross slip system. In this way cross slip arises 
although the dislocations never change their slip planes. 

This does not mean that cross slip does not also come about through 
dislocations moving from one slip plane into another—indeed this process 
seems to have been directly observed by Hirsch et al. (1956). Lf significant 
numbers of dislocations do change their slip planes, then, in a metal with 
short range order, there is a high probability that several dislocations 
will follow the same path, a process which will lead to clear cross slip. 
However, in the case of pure metals no noticeable difference exists between 
slip planes on which dislocations have moved already and those which are 
still undisturbed (except possibly after severe straining and near surfaces). 
Therefore in them dislocations changing their planes will do so more or 
less at random, and no sharp cross slip traces can be expected as a result 
of this process. Sharp cross slip can be produced in pure metals only 
under unusual conditions, such as, for instance, in the strongly inhomo- 
geneous stress field near an etch pit, where the two-fold action of 
R-dislocations, discussed above, takes place. 

In agreement with this idea, the expected distinct difference between 
the diffuse cross slip linking slip bands, and sharp cross slip lines near 
etch pits has been often observed in aluminium. Several examples of 
this can be found in the micrographs published by Wilsdorf and 
Kuhimann-Wilsdorf (1952) and by Noggle and Koehler (1957). 


§ 5. Some Par?vicULAR ADVANTAGES OF THE PRoposED THEORY 


At the beginning of this paper five sets of experimental results have 
been listed which appear to be in disagreement with the concept of the 
emission of dislocations from a random dislocation network. We shall 
now see that part of these problems can be solved with the theory proposed 
in this paper. 

(i) In order to make a quantitative comparison between the number 
of dislocation sources apparently active during the formation of an 
elementary structure, and the number of sources which can be expected 
through the condensation of vacancies, use will be made of measurements 
made on super-purity aluminium by Wilsdorf and Kuhimann-Wilsdorf 
(1952). The measurements had been done with great care but were 
published with the warning that ‘“‘ not too much reliance must be put 
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on such measurements” because the lines were not clearly resolved, 
their average height being below the resolving power of the preparation 
techniques available. 

Using the formula A=(1/d)[exp (—U/kT)] derived in § 3, and assuming 
that on the average the stacking faults are hexagonal with sides L we 
find n, the number of stacking faults per cubic centimetre of material, 
as n=(2-6L2d) Jexp (—U/kT)]. Each of the hexagonal R-dislocations 
represents four possible dislocation sources, so that on each of the twelve 
glide systems n’=n/3 sources per cm? can be expected. Taking, for 
example, a resolved shear stress of 1 kg/mm?, corresponding to a source 
length of about 10-4cm, and using U=0-76 ev, due to Bradshaw and 
Pearson (1957), which is quoted by Cottrell (1957) as the latest measurement 
of U in aluminium, we obtain n’=4-5 x 10!%/cem3. This value must be 
compared with the following measurements which represent the mean 
obtained from two different crystals stressed to 1 kg/mm*: mean distance 
between the elementary lines s=400 A, average length of lines D=17-5 p. 
The density of active sources can be found as 


1/ V =(4t7D?s)- 1 G-7 p< 1019; em?. 


In view of the several simplifications and uncertainties involved this 
agreement is most satisfactory. 

(ii), (iii) No complete answer to these points can be given in the present 
paper. It is believed that the strong dependence of surface structure and 
plastic properties on heat treatment often is based on the same phenomenon 
that was discussed above in connection with the «-brass structure, i.e. the 
changes which may occur on a slip plane as a result of the passage of the 
first dislocation or the first few dislocations, to move along that plane. 
It would appear that in many cases the changes in surface structure and 
plastic behaviour are not primarily due to any variations in the source 
density or their properties. A theory on this basis was outlined in the 
lecture mentioned before (Kuhlmann-Wilsdorf et al. 1957,) and it is hoped 
to deal with it more fully soon in a separate paper. 

(iv) All possible types of R-dislocations represent dislocation sources 
of the non-homogeneous type which cannot give rise to ‘ homogeneous ’ 
slip. This follows from a principle first given by Seitz (1952) that the 
area circumscribed by the projection in slip direction of a general 
dislocation ring—of which R-dislocations are examples—cannot be changed 
through slip of the dislocation ring or any parts of it. 

(v) The most important argument in favour of the present theory 
would seem to be the results obtained by Gilman and Johnston. Not 
only do they represent the most serious challenge to the theory of the 
three-dimensional dislocation network, but moreover they seem to be in 
excellent agreement with the generation of dislocations from sheets of 
condensed vacancies—stacking faults in the case of face-centred cubic 
metals. In fact Gilman and Johnston (1957) give as one of three possible 
explanations for their results that “small unseen dislocation loops might 
be present which grow to form Frank—Read sources when a stress is applied. 
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Although small loops on the glide planes would be expected to be unstable 
and hence not present, prismatic loops, which could not disappear by 
gliding might form by vacancy condensation (J.C. Fisher and F. C. Frank 
private communication) ’’. However, the authors then proceed to discount 
this possibility for the reasons that such prismatic dislocations ought 
to be eliminated through slow cooling and that they should be easily 
seen. 

While the first of these two reasons appears to be unsound the second 
is more difficult to judge by anyone not very familiar with the experimental 
technique used. It should be noted, however, that the necessary 
k-dislocations would have dimensions in the order of 1 micron only, and 
that they would hardly be distinguishable from the dislocation loops 
grown from them. Hence both reasons are probably not convincing 
enough to discredit the idea of R-dislocations as the origin of dislocations. 


§ 6. THE CHOICE OF GLIDE SYSTEMS IN VARIOUS CRYSTAL LATTICES 


Except for the satisfactory explanation of the fact that the active slip 
direction is always the most densely packed lattice direction, advanced by 
Nabarro (1952), no general and detailed theory on the choice of slip systems 
has so far been developed. In the idea of dislocation formation through 
condensed sheets of vacancies, presented here, the glide systems operating 
in any type of lattice are determined through the slip direction of the 
R-disiocations and through the orientations which the axes of the 
R-dislocations will tend to assume. If these parameters can be predicted 
then a theory of slip systems is obtained. 

At present our knowledge of vacancy condensation is scant, and no 
experimental data are available from which we could definitely deduce 
on what crystallographic planes it takes place, whether condensation layers 
are single or multiple, and what are the intermediate steps which lead to 
the nucleation of a fault. However, primarily vacancies are tiny holes 
in the material, and, with this idea in mind, we may make the following 
assumption. Vacancies, which accidentally meet during their thermal 
motion, form loosely bound groups of a few—say up to seven. The energy 
of such groups will be lowest if they take the shape of discs parallel to the 
most densely packed planes. If in the course of their thermally activated 
motion two or more such groups meet on the same plane they combine to 
form a dise-shaped hole large enough to allow its collapse, thereby creating 
an R-dislocation. Or it might even be that already the largest of the 
primary groups collapse. 

Hence we assume vacancies to condense preferentially on the most 
densely packed atomic planes, and the faults thus created to be bounded 
by the most densely packed directions in that plane. These, then, will be 
the directions of the axes of the R-dislocations, and, at the same time, will 
be the lines of intersection between the slip planes and the plane of the 
fault. The slip direction is found by considering what type of fault will 
arise when the two sides of a sheet of vacancies collapse, and what movement 
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will restore the correct lattice order. The imaginary path which any 
atom of one side of the sheet of vacancies will traverse with respect to the 
other side of the sheet, from the start of the collapse to the re-establishment 
of order, is the slip vector. In this way both, slip plane and slip vector, 
can be predicted, but it should be added that the slip planes thus found 
are not necessarily those which offer least resistance to the movement of 
dislocations. Therefore the predicted slip planes could be supplemented 
by planes into which the dislocations move through cross slip motion. 
Good reviews of glide systems found experimentally in a number of 
different lattices are contained in the book by Schmid and Boas (1950) 
and in a more recent review article by Maddin and Chen (1954). 

(i) The face-centred cubic lattice has been partly examined already. 
The vacancies condense in {111} planes forming stacking faults bounded 
by (110) directions. The movement due to the collapse of the sheets of 
vacancies, or in other words, the Burgers vector of the resulting sessile 
dislocation is 4<111). The movement which restores order, i.e. the 
Burgers vector of the half-dislocation, is 4<112), and the slip vector is 
the sum of the two movements, namely 4$[111]+23{112]=3[110]. 

Slip planes are all planes which contain the slip direction of the 
R-dislocation and a ¢110) direction lying in the plane of the fault. We 
saw (cee figure) that two of the three slip planes defined in this way are 
{111} planes. The third plane, so far neglected, is a {100} plane, a slip 
plane which is reported on aluminium at elevated temperatures. 

(ii) The geometry in the hexagonal lattice is rather complicated, and it 
ought to be investigated more closely than is done here. It is unique in 
that removal of one {0001} plane yields a stacking fault which cannot 
be removed by slip (Read 1953), and the same holds true for the removal 
of one {1010} plane. Consequently no slip dislocation can be formed by 
the amidenention of single layers of vacancies on either the basal plane 
or the prism face. The Burgers vector which results in both cases is 
(1013). This vector lies in the {1122} plane, recently identified a as 
slip plane in zinc by Bell and Cahn (1957), and includes a small angle only 
with (1123), the slip vector found by the same scientists. However it 
is extremely improbable that a dislocation with (1013) as Burgers vector 
could move because of the severe lattice disturbance which would thereby 
be created. 

The slip vector (1120), commonly observed as slip direction in all 
hexagonal metals, could arise through vacancy condensation on either 
{1120} or {1122}. The most probable slip planes associated with these 
planes are {0001}, the basal plane, and {1010}, the prism plane. 

(iii) In the NaCl structure, which is also the structure of lithium fluoride. 
vacancies would be expected to condense in {100} planes, and the resultant 
faults to be bounded by <100) directions. The removal of one {100} 
plane must be accompanied by an offset of ${100) because otherwise 
like ions would face each other. The slip vector thus becomes 


${100]4-$[010]—4/ 110], 


: tm) = hen 
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in agreement with experimental fact. The planes which intersect { 100} 
planes in (100) directions and contain the (110) slip direction are 110} 
and {100}, representing possible slip planes. {100} is the slip pint 
actually observed in this lattice type. 

(iv) The body-centred cubic lattice, which is geometrically identical 
with the caesium chloride structure, was recently investigated by Rachinger 
and Cottrell (1956) because of its interesting property that the slip direction 
found experimentally is (111) when the atomic bonds in it have metallic 
character, but is (100) is they are ionic. Rachinger and Cottrell succeeded 
in interpreting this dependence of slip direction in terms of partial 
dislocations. Another simple explanation follows from the present theory. 

The most densely packed lattice plane is {110}. Removing one {110} 
plane leaves a fault which will be rectified by an offset of either $¢100) 
or $<{110). The former of the two possibilities represents the smaller 
lateral motion and, therefore, must be strongly favoured in the body-centred 
cubic lattice. In ionic crystals with the same structure, on the other hand, 
an offset of $<(100) would cause like ions to face each other, which 
obviously is energetically impossible, and hence, the offset in that lattice 
will be $110). The resultant slip vectors become $[110]+4[001]=4[111] 
for the case of metallic bonds, and ${110]+43[110]=[010] for ionic bonds, 
in agreement with experimental evidence. 

As before, the slip planes are found as the planes which are common 
to the slip vector and the most densely packed lattice directions in the 
plane of the fault. In this case the (111) directions are the most closely 
packed, and two of these lie in any {110} plane, which, as we saw, is the 
plane of the R-dislocations in the body-centred cubic lattice as well as 
in the caesium chloride structure. In metals with (111) as slip direction 
both <111) directions in the plane of the fault define {110} as the slip 
plane, and this is in fact the dominant slip plane found experimentally. 
Also if <100) should serve as bounding direction {110} results as slip 
plane, while with (101) as boundary—the third densest lattice direction 
in body-centred cubic lattices—{121} becomes a possible slip plane. 
The third slip plane reported for body-centred cubic metals, {123}, cannot 
be explained in an obvious way. 

The failure to account for {123} as shp plane, and the rather remote 
probability that {121} would result through vacancy condensation might 
lend support to the contention of Maddin and Chen (1954) that the wavy 
slip lines of body-centred cubic metals, usually attributed to the simul- 
taneous action of the three types of slip planes mentioned, could be the 
result of intimate cross slip on {110} planes only. [It is also possible, 
though not very probable, that in the body-centred cubic lattice a plane 
other than {110} offers the least resistance to the motion of dislocations, 
and that, therefore, many dislocations change to that slip plane by means 
of cross slip—{110! being the only plane on which dislocations can be 
formed. 

For ionic crystals with (100) as slip direction again {110} is the slip 
plane defined by both (111) directions in the plane of any K-dislocations, 
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The two next densest directions, (100) and (110), would yield {100} as 
possible slip plane. Experimentally {110} is found as slip plane in the 
caesium chloride structure. 

The above theory of slip systems rests on the assumption that vacancies 
condense in single layers parallel to the most densely packed planes. If, 
on the other hand, we were to assume that the vacancies do not condense 
in this way, but such that they yield R-dislocations with the smallest 
possible energy the resultant agreement between theory and experimental 
fact would be much poorer. In’ all cases double layers normal to the 
shortest lattice vector give the smallest possible energy, so that the slip 
directions would be accounted for. However, probable slip planes are 
found in the order given: body-centred cubic... {112} and {101} ; 
MsCl... {100} and {110}; NaCl... {110} and {100}; hexagonal 
lattice... {0001}. Thus with the exception of the hexagonal lattice, 
the dominant slip planes always appear as second choice. 
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ABSTRACT 


Explicit formulae are found for the total drag on a class of obstacles placed 
in a uniform stream of an inviscid fluid when the stream function is given 
in terms of the velocity variables. The formulae are given for incompres- 
sible flow and for the exact equation of compressible flow. 


THERE are a number of problems in fluid dynamics for which an exact 
solution for the stream function Y can only be found in terms of the 
velocity variables. Such problems are said to be solved in the hodograph 
plane. It is one of the disadvantages of the hodograph method that 
the transition to the physical plane can seldom be effected in a simple 
manner. The object of this note is to show that the drag on any straight 
segment of a body can be found directly from the solution in terms of the 
velocity variables. 

The use of the hodograph method in incompressible flow dates from 
the time of Helmholtz and Kirchhoff. More recently it has been 
extensively used in the study of compressible flow, particularly in the 
transonie region where the usual linearization procedures break down. 
In incompressible flow Y is a harmonic function of w and v, the cartesian 
components of velocity. It is usually more convenient to work with 
the polar coordinates q and 6. For compressible flow we introduce a 
variable 7 defined by t=q?/q,,? where q,, is the maximum velocity attain- 
able by the gas. YW now satisfies a second order linear partial differential 
equation in 7 and @ and the solution of this obtained by separating the 
variables is of the form %,(7) sin (n6+-e,). (7) is a certain hyper- 
geometric function multiplied by a power of +r. The details of this 
procedure and the main properties of %,(7) are to be found in Modern 
Developments in Fluid Dynamics (1953). For our purposes we regard 
s,,(7) sin (nO+-e,,) as the compressible analogue of the harmonic function 
q" sin (rf--e,) which is the solution appropriate to incompressible flow. 

We are concerned with a class of problems in which an obstacle is 
placed in a stream which is uniform at infinity with velocity qo. @ is 
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measured from the direction of this uniform stream. The obstacle is 
such that its profile consists of segments along which either the magnitude 
or the direction of the velocity vector is constant. This may appear 
to be a rather severe restriction but in fact it is only for obstacles of this 
type that the boundary value problem can be exactly formulated and 
hence solved in the hodograph plane. We consider the drag on a straight 
segment inclined at an angle « to the free stream. For incompressible 
flow the velocities at the two ends will be denoted by qg, and q,. The 
stream function will be of the form 


P=)a,q" sin (nb-+e,). 


If p and p are respectively the pressure and the density, then the 
excess pressure over the pressure at infinity pp» is given by Bernoulli’s 
theorem in the form 


P—Po=2P%o —3PT- 
The total drag in the direction of the main stream is 


sin « | pds=L sin «(po+4pqo")—Fp sin « \¢ ds, 


where the integration is taken along the body contour and L is the total 
length. Now on the boundary Y is a constant which we can take to be 
zero and g=dd¢/ds where ¢ is the velocity potential. Hence 


i r (2 0d 
2 ds= |q d6= = dq. 
|@ ds=|q dg ) Taq 
Corresponding to Y=q” sin (nO+e,) we have d=—q” cos (nf+e,). It 
follows that the total drag is given by the expression 


MGinp SIN 4 COS (MOAPER) (9 n4t__g n+1) 
2(n+1) ; aa 


In this formula cos (nx+e,,)=+1 since sin (na+e,)=0. 

For a similar straight profile in compressible flow we denote by 7, 
and 7, the values of 7 at the two ends. The appropriate form of 
Bernoulli's equation is now 


L sin «(pp +4p402)-+ > 


p=p (1—7)P tt 


where p’ is the stagnation pressure, B=1/(y—1) and y is the adiabatic 
index of the gas. The velocity potential ¢ corresponding to a solution 
P=u,,(7) sin (n?+e,) 1s 


d=— = (1—r)-Fib,,'(7) cos (nO+e,) 


where the dash denotes differentiation with respect to 7. The contri- 
bution to the total drag arising from this term is found to be 


wenn fe fee H Soa aS ide 
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It is rather surprising that the integral in this term can be evaluated 
explicitly. This is achieved by substituting for %,,(7) in the original 
partial differential equation and integrating by parts. Corresponding 
to the solution 


Y= >a,,(7) sin (nO+e,,) 


which is to vanish on 6=« we find, after a certain amount of algebra, 
that the total component of the drag in the direction of the main stream is 


p'na,, sin « cos (nx-+e,,) 
Gil) 

The summation can be over any sequence of positive numbers » which 
does not include the value n=1. This restriction is of no significance 
as it can be shown that the presence of a term corresponding to n=1 in 
the solution requires the fluid to be at rest at infinity. Such a term 
cannot therefore arise in problems of flow past a body such as we are 
considering. 

A similar result may be shown to hold for Tricomi’s equation which 
is an approximation to the full hodograph equation and valid for small 
perturbations about a uniform flow at a Mach number of one. 


[272/2y,'(r)(1—7) +7 Mp, (7) {1+ (28-+ Ir}. 
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ABSTRACT 


Energy spectra and angular distributions of protons emitted from the 
isotopes *4Fe and °*Fe on bombardment with neutrons with an energy of 
13-5 Mev have been studied using photographic emulsions. The angular 
distributions of protons with energy greater than 7 Mev show a marked 
forward peak suggesting the occurrence of direct interaction. 


§ 1. INTRODUCTION 


Srupies of the inelastic scattering of protons with energies in the region 
10-30 Mev (Gugelot 1954, Eisberg and Igo 1954, and Cohen 1957) have 
indicated the occurrence of interaction processes other than the immediate 
formation of acompound nucleus. Brown et al. (1957), found evidence that 
direct interaction also occurs in (n, p) reactions at 13-2 Mev. To obtain 
more information about the direct interaction mechanism, similar studies 
of (n, p) reactions were undertaken by the present authors with improved 
energy resolution. 

The isotopes °4Fe and °*Fe were chosen for investigation because the 
effects of the Coulomb barrier should be similar, whereas the (n, p) cross 
sections should differ considerably in magnitude due to the differences in 
the @ values for the (n,p) reaction. The relative contributions to 
proton emission of compound nucleus formation and of the direct inter- 
action might be expected to be different for each isotope (e.g. Brown and 
Muirhead 1957). 


§ 2. EXPERIMENTAL PROCEDURE 
2.1. Exposure and Examination of the Photographic Plates 


The separated isotopes consisted of layers 1 cm? in area and 7 mg/cm? 
in thickness electroplated on to gold foils. The °*Fe separated isotope 
contained 99-9°% 5*Fe, while the *4Fe separated isotope contained 67%, 
54Me and 32% °**Fe. The separated isotopes were placed in contact 
with the surface of an Ilford C2, 4in.x1lin.x400,, nuclear research 
emulsion and covered by a gold strip of size 4in.x1in.x400py. The 
gold, which has a very small (n, p) cross section, was used to shield the 
surface of the emulsion from external sources of protons. 

Six emulsions were exposed simultaneously mounted hexagonally 
around the neutron source using the geometrical arrangement described 
by Brown et al. (loc. cit.). Two of the emulsions were in contact with 
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54Fe targets, two with **Fe, and two were exposed without an iron 
isotope between the emulsion surface and the gold strip. 

The neutrons, produced by the T(d, n)*He reaction, were incident on 
the isotopes at an angle of 30°--2° to the emulsion surface and had an 
energy of 13-5-L0-1 Mev. 

The emulsions were scanned using X50 oil immersion objectives and 
x10 eyepieces for tracks entering the surface of the emulsion. All 
singly charged tracks were assumed to be protons. (Deuteron emission 
was expected to be small as a result of unfavourable @ values.) Only 
tracks whose angle of dip in the emulsion was between 10° and 45° were 
accepted for measurement. These limits were taken to minimize the 
possible energy loss in the isotope layer and to ensure that the majority 
of the high energy protons stopped within the emulsion. Hence, protons 
making angles greater than 140° with the direction of the incident neutrons 
were not accepted. 

The range, angle of dip and the angle in the plane of the emulsion 
with respect to the direction of the incident neutrons were measured for 
every track accepted. The energies of the protons were obtained using 
the range energy tables of Rotblat (1951). 


2.2. Subtraction of Background 


The majority of the background protons are due to the interaction of 
neutrons with the emulsion nuclei within a distance of 1-2 microns from 
the emulsion surface. The tracks produced by these protons cannot be 
optically resolved from those protons entering from outside the emulsion. 
‘The only other source of background is the gold shield. 


Table 1 


Number of Number of 
Isotope protons Background | protons from 
observed isotope 


54Be 1823 1502 


6He 2479 


The emulsions covered directly by the gold strips were used as control 
plates ; the distributions in angle and energy of the protons observed 
in these plates were subtracted directly from those found in equal areas 
of the emulsions in contact with the separated isotopes. 

This procedure neglects the effect on the range of the protons from the 
gold strips of their passage through the isotope layer. However, since 
these protons comprise only about 20°% of the total background and the 
energy loss in the isotope layer is small for many of the events, this 
procedure is considered justified. 

The total number of protons observed entering the emulsion surface 
for each isotope and the background subtractions are given in table 1. 
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2.3. Correction for Energy Loss 

The protons originating within the separated isotopes lose energy in 
their passage through the isotope layer to the emulsion surface. Allowance 
was made for this loss by adding to the observed energy of each proton the 
energy it would have lost in traversing half the thickness of the isotope 
layer, taking account of its observed angle of dip with respect to the 
emulsion surface. Thus, a proton having a dip of 30° and having a range 
in the emulsion corresponding to an energy of 4 Mev would be considered 
to have an initial energy of 4:3-.0-3 Mev. 


2.4. Determination of the Neutron Flux 


The neutron flux was obtained by scanning a known volume of emulsion 
for protons resulting from the collision of the incident neutrons with the 
hydrogen present in the emulsion. These recoil protons, which produce 
the majority of the tracks found within the emulsion, are easily identified 
by applying the relationship, 


Ep=(13-5+0-1) cos? @ 
where Hp is the energy of the recoil proton and @ is the spatial angle 
between the direction of the recoil proton and the incident neutron. 


§ 3. RESULTS 

The total cross section for the emission of protons from each isotope 
was obtained by extrapolating the angular distribution to 180°. The 
values obtained were 600--60 mb for °4Fe and 95+20 mb for ®**Fe. 

The angular distributions of the protons emitted from °*He and °°Fe for 
three different ranges of proton energy are shown in figs. | and 2. 

The energy spectra of the protons emitted from °4Fe and °*Fe are 
shown in figs. 3 and 4. The three distributions given in each figure 
cover the angular ranges 0°-60°, 60°--90°, and 90°-140°. 

The solid line curves in figs. 3 and 4 are the proton distributions 
expected for the (n, p) reaction on the basis of the statistical theory of 
Weisskopf and Ewing (1940). In these calculations the Coulomb barrier 
penetrabilities given by Blatt and Weisskopf (1952) for a nuclear radius 
of 1-3 10-13 Al3 cm were used together with the formula given by Lang 
and Le Couteur (1954) for the density of nuclear levels as a function of 
excitation energy. The calculated energy distribution for °*Fe was 
normalized to the observed distribution for the angular range 90°-140°. 
The curves shown for the other two angular ranges were obtained by 
assuming that the angular distribution of protons due to this process is 
isotropic. For 54Fe the Coulomb barrier should be almost identical, but 
the number of protons emitted with low energy is considerably greater 
than for 5*Fe; this is presumably due to the (n, np) reaction (see § 4). 
Thus for *4Fe the calculated (n, p) distribution was normalized to the 
observed distribution for proton energies greater than 4 Mev in the 
angular range 90°-140°. The broken line was drawn through the observed 
distribution for protons with energies below 4 Mev in this angular range. 
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The curves drawn on the other two angular intervals were obtained by 
assuming that protons emitted by these processes are emitted isotropically. 


§ 4. Discussion 

The angular distributions of the protons emitted from both °*Fe and 
56Ke are approximately flat at angles greater than 90° for all energy 
intervals. To illustrate any departure from isotropy at forward angles, 
the dashed line was drawn through the mean of the experimental points 
at backward angles. The protons emitted from both isotopes with 
energies less than 4 Mey are approximately isotropic. For proton energies 
between 4 and 7 Mev there is a tendency for the distribution to peak at 
forward angles and for proton energies greater than 7 Mev there is a 
marked peak at angles below 60° for both isotopes. 

The preferential emission of protons in the forward direction at high 
energies cannot be explained by the statistical theory for the decay of the 
compound nucleus but would be expected if protons were emitted as an 
immediate result of the collision of the incident neutron with a proton 
of the nucleus. 

The large number of protons with energies below 4 Mev emitted by 
54Re is probably due to the (n, np) reaction which is favoured relative 
to the (n, 2n) reaction due to the differences in @ values for the two. 
reactions as shown in table 2. 


Table 2 
: : Maximum energy of 
Isotope SEAS ORY emitted proton 
(n, p) | (a, np) | (n, 2m) | (n, np) (n, p) 
54Re | +0-34] — 9-2] —13-8 4-3 13-9 
56He | —2-93 | —10-2 | —11-2 3:3 10-6 


Thus in °*Fe the (n, 2n) reaction cannot occur when the energy of the 
incident neutron is 13-5 Mev, while the (n, np) reaction is possible for 
proton energies from 0 to 4-3 Mev and so has only to compete with the 
(n, ny) reaction. 

For °°Fe the (n, np) reaction is only favoured over the energy range 
0-1-0 Mev where proton emission is prevented by the Coulomb barrier. 

Figures | and 2 show that protons emitted from °4Fe with energy 
below 4 Mev have an isotropic angular distribution which is consistent. 
with an (n,np) reaction involving the decay of a compound nucleus. 
The good fit of the higher energy protons to the calculated energy spectrum 
for the angular range 90°-140° for both isotopes suggests that these 
protons also result from the decay of a compound nucleus. 
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For both °*Fe and °*Fe there are more high energy protons at forward 
angles than would be expected from compound nucleus decay. These 
protons are probably due to a direct interaction process. Their energy 
spectrum is approximately flat from 4 Mev to the maximum energy allowed 
by the Q@ value for the (n, p) process. The statistical accuracy of the 
results are such that the presence or absence of structure in the energy 
spectrum, similar to that found by Gugelot (1954) and Cohen (1957) in 
the inelastic scattering of protons, cannot be confirmed. 

If the isotropic part of the angular distributions are assumed to be due 
to compound nucleus processes and the forward peaking to a direct 
interaction, we get the cross sections for these processes given in table 3. 


Table 3 


Cross section in millibarns 


Isotope 
Compound nucleus Direct Total Allan 
(n, p) (n, np) | interaction 8 total 
*4Fe 306 224 70 600+-60} 680 
, 95+20| 190 


The direct interaction cross section obtained in this way does not 
include any isotropic component and so may be an underestimate. 

The reaction denoted (n, p) in this table includes both the reactions 
(n, py) and (n, pn). The reaction (n, pn) does not suffer as a result of 
competition with the (n, 2n) reaction, and is possible for proton energies 
in the range 0-4-3 mev for °4Fe and 0-3-3 mev for °*Fe. In these energy 
intervals it will dominate over the reaction (n, py) but at higher energies 
the protons arise entirely from (n, py) reactions. Hence a comparison 
of the cross section for proton emission above these energies with the 
(n, py) cross sections determined by activation method can be made. 
Paul and Clark (1953) obtained a value of 97-12 mb and Forbes (1952) 
a value of 124--12 mb for the reaction (n, py) in °*Fe at 14mev. The 
present results give a value of 934-20 mb for the (n, py) reaction at 
13-5 Mev. 

Allan (1957) has measured the energy spectrum of protons emitted 
from °4Fe and ®*Fe at an angle of 34°-+-20° to the incident neutron direc- 
tion. The total cross section he obtained is shown in the final column 
of table 3. Allan obtained these figures by assuming an isotropic angular 
distribution, and thus probably over-estimated these cross sections by 
about 30°. 

Austern et al. (1953) suggested that proton emission was enhanced by 
interactions at the nuclear surface. The shape of the angular distri- 
bution was characterized by the allowed values of orbital angular 
momentum with which the neutron can be captured by the nucleus. 
This theory does not, however, give absolute cross sections. 
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Brown and Muirhead (loc. cit.) consider a model which assumes that 
nucleon-nucleon interactions can occur throughout the nuclear volume 
and eject protons directly. They calculate the cross section for the 
(n, p) direct interaction at 14 Mev to be 55 mb for **Fe and 30 mb for 5¢Nes 
These values are in reasonable agreement with those obtained here. 
Their energy spectra and angular distributions, while showing the main 
tendencies of the present results, are not in detailed agreement. 

The surface interaction calculations of Elton and Gomes (1957) are 
also in qualitative agreement with our results. 

The present results do not determine the details of the direct inter- 
action process sufficiently well to make significant comparisons with 
these theories. 


Fig. 5 
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Proton spectra of 54Fe and 5*Fe corrected for Coulomb barrier penetrabilities. 


Recently Brown and De Dominicis (1957) have formulated the problem 
of inelastic scattering of nucleons by nuclei in terms of nuclear dispersion 
theory. They show that expressions for the contribution of both the 
‘direct interaction’ and of the ‘compound nucleus’ to the inelastic 
scattermg may be obtained by this method. 


In fig. 5 
N(E 
log, ( en) ) 
. Oo,, Hip 
is plotted against Hp, where N(Hp) is the number of protons observed 
with an energy #p in the angular range 90°-140°, and c,,, is the cross 
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section for the formation of a compound nucleus as given by Blatt and 
Weisskopf (loc. cit.) for a nucleus of radius 1-3 10-13 A” em. The 
experimental points are consistent with a straight line for 5*Fe and two 
straight lines for *4Fe. On such a plot a straight line indicates that the 
density of nuclear levels varies as exp(#) where H is the excitation energy 
of the residual nucleus, and its slope is equal to 1/7’ where 7’ is the nuclear 
temperature at the average excitation energy. The figure illustrates the 
occurrence of the (n, np) reaction in °4Fe at energies below 4 Mev and 
shows that a temperature of 1-2 Mev gives the best fit for the (n, p) reaction 
in both isotopes. 


§ 5. CONCLUSIONS 


The major contribution to the (n, p) reaction at 13-5 mev in both 54Fe 
and °®Fe is that made by the decay of the compound nucleus. In both 
isotopes considerable direct interaction effects have been observed. The 
protons resulting from direct interactions are emitted preferentially at 
forward angles and have an approximately flat energy spectrum. 
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ABSTRACT 


The motion in the boundary layer when a body of revolution is started 
to move along its axis of symmetry and simultaneously to rotate about it 
with a velocity varying as the nth power of time, is investigated. The 
corresponding problem for an impulsive start was tackled by Illingworth. 
It is found that the separation is delayed by acceleration. Also, the time 
that elapses before separation first takes place is lessened by increasing the 
velocity of spin. 


§ 1. INTRODUCTION 


THE aim of the researches on unsteady boundary layers has been to- 
determine the time that elapses after the start when separation first takes. 
place. The method of Goldstein and Rosenhead (1936) has proved quite 
satisfactory to this end, and describes the early stages of the development. 
of the boundary layer almost accurately. The method has been extended 
by Illingworth (1954) to consider the motion due to an impulsively 
started spinning body of revolution. Squire (1954) has discussed the 
growth of three-dimensional boundary layers on the same lines—again a. 
case of impulsive start. Wundt (1955) has studied the early stages of the 
boundary layer growth on yawed cylinders when the potential velocity 
varies as the nth power of time. The case of impulsive start can be 
obtained as a special case. It is worthwhile considering the Illingworth’s. 
problem, viz., the motion due to a spinning body, the translatory and 
rotatory motions being acquired from rest by an acceleration varying as. 
some power of time. The aim is to show that it takes more time for 
separation to occur in the accelerated case than in the case of impulsive: 
start. The extension of the present analysis to Squire’s paper is straight- 
forward. 


§ 2. EQUATIONS OF MorTION 
As in Illingworth (1954) we choose the independent variables x, y, and t 
such that ¢ denotes the time from the start of the motion, 2 measures the 
distance from the nose of the body along any meridian curve, and y is the 
distance along the normal to this curve. This means that the frame of” 
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reference moves forward with the body but does not rotate with it. If 
U, v, w, are the velocity components in the directions of x, y and the 
‘transverse direction respectively and r(«) is the radius of the body at the 
station x, the equations of motion are 


0 0 
ay lat ios (rv) =0, re ie eee ty aL) 
du Ou du wdr aU OU oy . 
Gr ncyo hide 4ae. eer aye. 
Ow Ow Ow . uwadr O2w 


where U(x, t) is the velocity of the main stream. It is supposed that the 
body and the fluid are at rest at <0. At t=0, the body is given a trans- 
datory velocity and a rotatory velocity varying as ¢” so that the boundary 
conditions to be satisfied are 


v——w—U" for ¢<0vand foriall-y> §. .°. .'. (4) 
u=v=0, w=2rt" at y=0 

== Uf = f ’ eae 
wet) =P (e)).'t nae + for t20 (5) 
w=0 


We use the Stokes’ stream function % such that 
reg Ors ae 1 ob 


T Oy’ 7 Ox” 
and choose 
n=y/2/(vb), Lt tA ih A Cee at Be Set, ae Gantt ats (7) 
= 24/(rt) y pee teal) Tee, 4); oe sk le (2) 
0) 
w= > trt4(n-+1) OT (x, 7). gl ae te 4 Pe en |) 


Substituting in equations (2) and (3), equating coefficients of like powers of 
#, and rearranging we obtain linear partial differential equations to deter- 
mine Y,,,and ®, ,. These can be reduced to linear ordinary differential 
equations if we choose 


Y,, (% n=r Pin i | 
ae (%, =r Lé- Forl ie ~ fool ) + 7Q? & fasln |, 


r a 
P2aP dr 
%,,abesm—r | P* TS fal) + a fae eg Fal 


) 
P32 G2 P3 /dr\2 dP dr 
a =) af os(n) + we (z fas(n)+ 2?r dx di 198™) 
() 


dr 


+:2°Pr = wed 37 


) 
+287 (2) flo) | 
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dP 
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a°P 


®,, ern) =| PT gaan) +r 


2 
1p a 
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P2 
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dP\?2 
(Z) 9ao(n) +P atte Y33(7) 
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aa Joly ee Kp ” gaol | 


dP dr 
(11) 


2 dr 
Te) dasln)+ 2? 5 ea) 


If the differentiations with respect to 7 be denoted by dashes the resulting 


equations become 


Si" +2nf,"—4nfy'=—4n, x» tit aateeeelie 
fas” +2 fos” —4(2n+ 1) fa: = 48255 - , (126) 
fas” +2 fos" —4(82+2)fa::=483 - - - + - (12 c) 
9x +2ng1'—4ng1=9, (13 a) 
Jos +29 oi —4(20+ 1) o;=4ta:, (13 6). 
Yai +293; —4(32+ 2)93,=4t3;, (13 ¢) 


where 


Safi?’ —1, 

Sse=—fifr” 

S9g=—91 

Ss3i=Siter —hi'fer 

S832= 2h for’ —Sfifer’ hfe, 

S33=3fi foo 
— 2f1" feos 

Ssa=fi'foo’ —Si'foo: 

S35= —f1'foo' fife’ 


S36=S1 fos files’ — 2919 a1) 
S37=S1'fos —Si' fos, 
S33=Si'fos’ Silos’ —2f1'fos— 


hisses <ilai. 21 


9 
~919 29: 


toy= fifi 
too=2f1'91 Si 
ea —feigr's 
a= —fiJor fair’; 
too=2fy Jot for Gi thi G22—for9y 
—fi921' S19 22 — 2fa291'» 
tsa=f1'Jo2—f 291 > 
tss5=2fo0 Gi tsi J22—S19 20 > 
tsg= —fo391'; 
ts7=2fo3 J1—f e391 - 


(14) 


§ 3. SOLUTIONS OF THE EQUATIONS 
Equations (12) and (13) can be written in the form 


f" +2nf’—4mf=F (15), 
The two particular integrals of (15) are 
p> 2, Are 
m—1 
q=PE+ & Baran. BY, (16) 
r=0 
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where H=erfn=2x | exp (—£4) dé, and a=1/4/r 
0 


2(m—r) 


Bana GED) At 
2(m+r)Bo,1—2rA og, 
Boa Sein aelyo et bane 
Hence the complete solution of (15) is 
f=Ap+ Bq+p ju eal + lee d (17) 
= n, ee 
J p'q—py er: ==0'q 


where 4 and B are arbitrary constants to be determined from the boundary 
conditions 


filn=f1' ()=fisln)=fis'(n) | 
eis: Jij(n)=9. es at 7=0, cog tLS) 
: day lek cae H at 7= 0. (19) 


However, the values of the functions at the wall and hence the skin friction 
and the separation point cannot be obtained in a general form. We 
therefore consider the particular case of uniformly accelerated motion 
(n=1). The case of an impulsive start can also be obtained by putting 
n=0. 


§ 4. UNIFORMLY ACCELERATED MorIon 
For the case of uniformly accelerated motion n=1. The solutions of 
eqns. (12a) and (13 a) satisfying the conditions (18) and (19) for this 
particular case can be obtained as 
fi=n(1+ 397) E + 3a(1+77) exp (—7?)—F(a+7*), . ~ (20) 
91=— (1+ 2n?)H—2no exp (—n?)+1+2y%  . . . . (21) 
With these values of f, and g,, the values of s,; and ¢,, can be determined. 
The solutions of eqns (12 6) and (13 b) for this case (n=1) become 


1 1 1 2 1 16 3 2 
aa Faas lak a7 +H E migtlt 31 1) a5) 
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~ The functions s3; and t,, can now be determined, and eqns. (12) and 
(13 c) can then be solved. But because of the cumbersome algebra that 
sets in, the solution is given only upto the second term. 

The formulae for the skin friction, then, take the form 


0 
Ty (=) =2P(ax)a4/(upt) ite 0:42613 et -+-0-04008 i ih 
oy y=0 r dx 
r dr 
96.022 — — Ara eta wet cey l 
+0:21396.Q22 5, — a i (27) 
ew Bey ; dP 
7.—=—}t (Fa). = 2Qary/ (upt) J te oe = 
+0-21217 xs a +. aoe ee et bo (28) 


§ 5. SEPARATION POINT FOR A SPHERE 
We can now apply the results to bodies of given shape. In the case of a 
sphere we have 
Olax, t=. P(a)=t . 4A sin (z/a), 

Bir SUN a cme ee es ein en) oo eae (20) 
where 4 is the uniform acceleration with which the fluid moves relative to 
the sphere, and a is the radius of the sphere. Defining a non-dimensional 
time = Ai?/2a and B=a?.Q?/A? we have 
7 = 1-692574/(upt) . A sin (x/a) [1+ €(1-39863-+ 0-285288) cos (#/a)-+ . . .] 

(30) 

T -=1-128382/(upt) . a sin (a/a) [14-0-75675E-+ .. .]. fo GRY) 

Hence the separation first occurs at v=a7, i.e. at the rear stagnation 
point, and the value of € is given by 

1—€ (1-39863-+0-285288)=0. . . . . . (82) 


Since € denotes the distance through which the liquid particles at infinity 
have moved in time ¢, it can be compared with €, obtained by Illingworth. 


& 


0 0-71499 
(Del 0-70069 
1 0-59386 
10 0:23521 


Thus the table shows that €>€,, i.e., the liquid particles at infinity have 
gone much farther off in the case of accelerated start than in the case 


P.M, M 
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of impulsive start when separation first takes place. Also, as in the case 
of impulsive start, decreases as 8 increases showing that the separation 
occurs sooner the greater the velocity of spin. 
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Ty the zone-melting process of Pfann (1952), molten zones of constant 
length are caused to pass along a solid ingot of material, which it is desired 
to refine with respect to impurities whose distribution coefficients k between 
solid and liquid differ from unity. The coefficient & is defined as the ratio 
of impurity concentration in the solid to that in the liquid from which the 
solid is freezing. It is of interest to determine analytically the ultimate 
distribution of a given impurity in the ingot after the passage of an 
infinite number of molten zones. 

We make the usual assumptions of perfect mixing in the liquid and 
negligible diffusion of impurity in the solid; the value of the distribution 
coefficient / is assumed independent of concentration. Since the ultimate 
distribution y(x) is unaffected by the passage of a further molten zone, 
we may write down the equation (Pfann 1952) 


yla)=k[ y(€) dé, eT ate) 28 (1) 


where the distance x is measured in units of the zone length, and the zone 
moves in the direction Oz. 

Let the ingot be N zones in length, where N need not be integral. In 
the last zone-length of the ingot to solidify, i.e. for N—l<v<N, the 
length of the molten zone decreases progressively to zero ; accordingly for 
the last zone eqn. (1) is replaced by 


rN 
va=yof Od ...--- ®) 


—aX 4 


of which the solution is 
yl) = yo(u)=y(N—2)*1, N-1<e<N, . . . (3) 
where the constant y, is the impurity concentration at s=N—lI. . 
The distribution (3) is a boundary condition which has to be taken into 
account in the solution of (1). For example, the ultimate distribution 
Y,.(x) in the range N—2<x<N—1 is given by 


~N-1 ro+l 
yrale—k | aal€) ER] yorl) dé. 
x N=1 
From this equation we derive 
(D+k)yq2(%)=kyos(e+1) ee ee) 
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of which the solution is 
y4a()=A 12 exp (—kar)-+k exp (kx) | exp (leer) yor(e+1) de. 
With the aid of eqn. (3), integration successively by parts yields 
Y19(%)—= Ayo exp (—he)—y,(N—1—w)*M{1, 14h, b(N—1—2)}, 
No2<eeN—=1) eee 


a(a+1 
where M(x, 8 2)=14+ 4 co gen’ Pete 


is a confluent hypergeometric function, and the constant 
A .=Yy, exp {k(N—1)}.  ) jeu eas ae ee 
By a derivation in succession of the distributions Y1, Y19, Y23 - - - in the 
1, 2, 3,... zones from the end of the ingot, we find for the distribution in 


? 


the range N—(r+1)<x<N—,r, 
Yr, rt i(4)=A, p41 exp (—kx)+k exp (—kz) i exp (kx) y,_1, (a+1) da 
=A, .., exp (—hkx)-+-A,, , ke exp {—k@+1)}-+ -.. 
kerk 
(p= 1)! 
+khry, exp (—ke) | — | exp (kx)(N—r—ax)*-! dx, | 
where the constants 


Ay, or=ue xp {HN —)}— | Ay a, MVP) exp (—2) 


. +Ay— | («+r—2)"-l exp{—k(x+r—1)} (8) 


ke 
ries) 7-1 5] (Near) exp (—2k)+ 2S 


kr-1 
Cay ae aeeP {—@—np |. * ee 


and y, is the value of y(v) at ~=N—r. Term by term integration yields 
the result 


| es | exp (kax)(N—r—a)"1 dex7= (—1)* exp (kz) 


k(k+1)...(k+s—1) 

x (N—r—a)**""1 M{s, s+k, k(N—r—x)} (10) 
and this together with eqns. (8) and (9) gives the desired analytical 
expression for the ultimate distribution : 


k2 
Yr, rri(%)=exp tk(x—r) } | kya 51 (X—")Yr—2— 


(—k)'-1 Stas 
mi =n) vn | (11) 


—k)" 
+ EI ESD) A hb HoH} 
0<(y—r)<1_ r=1,2,...9| 


where we have put 
XN Sta as 5 ok ee 
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and v is the integer such that 0<(N—v)<1. It is readily seen that the 
ultimate distribution has a functional form which is independent of the 
length NV of the ingot. 

In the figure is shown the distribution (11) for the last four zones of an 
ingot, where the value k=0-1 has been chosen. This is in good agreement 
with the ultimate distribution calculated by Burris et al. (1955), who used 
a numerical method of solution of the ingegral equation. 


chia. Sha = “ao NA) N 


—> X 


The ultimate distribution of a zone-melted impurity, for which k=0-1. 
(————_) according to eqn. (11). 
(— — —) Pfann’s solution. 


By neglecting the effect on the ultimate distribution of the distribution 
in the last zone, Pfann (1952) wrote down as solution of eqn. (1) 


y(xz)=A exp (Az), 2 oh epee. Baa Gs) 
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where A is a constant, and \ is the non-zero root of the transcendental 
equation 

A=k (exp A—1). . i ee ee ee 
This distribution is also shown in the figure, with arbitrary constant A. 
The ultimate distribution (11) appears to approach Pfann’s solution as 
one moves further from the last zone; however, the author has been 
unable to demonstrate this analytically. 


Note added in proof.—The author is grateful to the referee for drawing 
to his attention a paper by I. Braun (1957, Brit. J. appl. Phys., 8, 457) 
on the ultimate distribution in the zone-melting process, published 
after the above was written. 
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ABSTRACT 


An investigation has been made of the long-lived ground state of 2°Al. 
The gamma ray spectrum has been determined with a scintillation spectro- 
meter, and the half-life calculated both by comparing the ratio of the masses 
of 76Al and ?7Al in the source and by comparison of reaction cross sections, 
the value obtained being 0-8 x 10° years (130%). 

The positron spectrum, measured using an overcompression cloud chamber, 
is consistent with that expected for a unique second forbidden transition, thus 
helping to confirm that the ground state of ?®Al is a 5+ level, decaying to a 
2+ excited level in **Mg. 

An Appendix describes the radio-chemical separation procedure for 
extraction of the active aluminium from the magnesium targets. 


§ 1. InTRODUCTION 


ALUMINIUM-26 is one of the odd—odd nuclei of mass number 4n-+-2 in 
which states of isotopic spin 7’=0 and 7'=1 are expected to lie close 
together, one forming the ground state. The indentification of these 
states depends on accurate determinations of Q-values for nuclear reac- 
tions and on a study of the cascade transitions in (p, y) capture reactions. 
Superallowed beta decay is expected between one or both of these two 
states and the J=0, 7'=0 state ground of the even-even isobars. 

In the case of mass 26, arguments such as these (summarized by Kluyver 
ct al. 1954 and by Kavanagh et al. 1955) suggest a decay scheme whose 
main features are shown in fig. 1. The two odd nucleons in ?®Al are 
expected to be in the d;,. shell and could couple to give, among others, 
states with (J, 7')=(0, 1) and (5, 0) of which the first would be the analogue 
of the 0* ground state of ?*Mg. ‘The 6-7 sec positron decay, first observed 
by Frisch (1934), is assumed to take place from a 7'=1 state identified 
in the (p, y) reaction at 230 kev excitation, the ground state decaying 
much more slowly to excited states of the stable isobar **Mg. The 
ground state decay was first reported by Simanton ¢é al. (1954), who 
detected positrons of energy about 1 Mey and gamma radiation of energy 
about 1-8 Mev, the half-life being estimated as about 10° years. 
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Handley and Lyon (1955) observed the spectrum of gamma radiation 
from 2°Al and Laubitz (1955) obtained evidence from the positron 
spectrum for the unique second forbidden nature of the ground state 
transition shown in fig. 1. The present work was undertaken with the 
object (a) of confirming the lifetime observed by Simanton et al, and 
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(6) of further establishing the energy spectrum of the positron decay. 
Since the work was finished the results of an extensive investigation of 
the lifetime have been published by Rightmire e¢ al. (1957) who established 
also a 16°, K-capture branch in the decay. Reports of decay scheme 
studies have also appeared from Ferguson (1957) and Johnson and 
Moffat (1957). 


§ 2. PREPARATION OF THE 26Al SOURCE AND ESTIMATE OF HALr-Lire 


Magnesium targets were bombarded by a 15 Mev deuteron beam in 
the Nuffield cyclotron of the University of Birmingham. 2°Al was 
produced by the reactions ?°Mg (d, n) and ?6Mg (d, 2n) and at the same 
time ?2Na was formed by the reaction 74*Mg (d, «). Both ?2Na and 26Al were 
extracted at the Radiochemical Centre, Amersham, the procedure for 
the 2°Al being described in the Appendix to this paper. The 2&Al source 
was placed close to a sodium iodide scintillation spectrometer and the 
gamma-ray spectrum was recorded, the main features, shown in fig. 2, 
agreeing with those found by Handley and Lyon. The half-life of the 
*°Al was estimated by calculating the efficiency of the scintillating crystal 
and by comparing the total gamma disintegration rate with the disinte- 
gration rate of the *Na made in the same bombardment. In another 
experiment the ?®Al activity was compared with the yield of 34Cl obtained 
by deuteron bombardment of a sulphur target. Estimates of lifetime 
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made by these methods depend on the assumption of equal cross sections 
for the reactions compared. The results obtained varied between 1:3 
and 4-6 x 10° years and were clearly inferior to those of Rightmire et al. 


(1957) who quote, as the mean of many observations, a half-life of 
1-1+0-1 x 10° years. 
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In order to avoid the necessity for making assumptions about the 
reaction cross sections a source of °Al containing a known weight of ?’Al 
carrier was analysed in a solid-source mass spectrometer. The ratio 
26 A1/27Al was found to be 2:45+0-37x10-°. From the gamma disinte- 
gration rate observed with the scintillation spectrometer the half-life 
was found to be 7',;,.=0-8+0-2x 10° years. This value, although less 
accurate, agrees with the results of Rightmire et a/., and confirms that no 
large error arises from the assumptions made in the ?Na/?°Al comparison. 


§ 3. THE Positron SPECTRUM 


The activity of the best available source was too low for the positron 
emission to be studied with a normal type of magnetic spectrometer and 
it was decided to use a 9 in. overcompression cloud chamber (Walker 
et al. 1956) fitted with a magnetic field of 500 gauss and filled with 
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hydrogen at atmospheric pressure. A suitable mounting for the source was 
made by spreading the active aluminium oxide as evenly as possible 
over a 5 micron aluminium foil, to which it was affixed by zapon. 

To check that the overcompression chamber was a suitable apparatus 
for such measurements, a source of Na, which is known to yield a 
Kurie plot of allowed shape, was mounted similarly to the aluminium 
and its positron spectrum investigated. Compton electrons arising from 
the accompanying gamma rays could be distinguished by their charge 
from the emitted positrons. No account was taken of the small number 
of positrons arising from pair production in the source, but this should 
not be a serious error, for such positrons would be of low energy, and 
would appear in that part of the spectrum already distorted due to 
source thickness. Within the statistical deviations, the Kurie plot was 
found to be linear from 150 kev to the end point of 540 kev, provided 
that only tracks lying within 25° of the horizontal were accepted for 
measurement, thus demonstrating that the overcompression chamber 
was suitable for investigating the activity of weak beta-emitters. 

With the ?°Al source in position it was found difficult to maintain 
satisfactory operation of the chamber, owing to the presence of an 
impurity which contaminated the filling gas. This meant that many of 
the expansions were of poor quality, and if tracks observed on all photo- 
graphs had been measured there would have been an excess of low energy 
positrons, because the lower the energy the greater the ionization and 
hence the greater the probability of seeing the track. A certain minimum 
quality of acceptance had thus to be defined, and in spite of this being 
largely a subjective criterion, two observers with independent criteria 
obtained similar positron distributions within statistical errors, although 
in one case 50°, more tracks were measured. There was rarely more 
than one measurable track per expansion, and less than 1000 suitable 
tracks were obtained from 4000 photographs. 

The positron distribution obtained is shown in fig. 3. The allowed 
Kurie plot (fig. 4), although not departing greatly from linearity, yields 
a maximum positron energy well in excess of the value of 1-16 (-- 0-02) Mev 
known from the reaction data. After application of the unique second 
forbidden shape factors (Konopinski 1955) the best straight line through 
the experimental points does yield the correct end point, although the 
statistics are not sufficiently good to establish this type of transition 
with certainty. Taken together with the evidence of other workers, the 
experiment supports the decay scheme given in fig. 1. 
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Extraction of ?®Al\ from Residues from Deuteron Bombarded 
Magnesium Targets 

The starting material was residues remaining after ?2Na had been 
extracted from deuteron bombarded magnesium. For extraction of the 
2Na the active surface had been dissolved off in acetic acid, and from this 
acetic acid solution the ?2Na (together with added carrier sodium) was 
precipitated by addition of excess uranyl acetate and ethanol. The 
residue therefore contained large quantities of magnesium, some uranium 
and its decay products, and traces of 22Na. 

1. The bulk of the ??Na remaining together with the bulk of the 
uranium was precipitated by addition of the correct quantity of sodium 
chloride solution. 

2. To the supernatant was added aluminium chloride solution 
(= 15 mg Al) and ferric chloride solution (= 100 mg Fe). The solution 
was evaporated to remove the bulk of the ethanol. 

3. Ammonia and ammonium chloride solution was added in presence 
of tannic acid. The precipitate contained the aluminium and _ iron, 
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together with a trace of uranium. The precipitate was washed, dried 
and ignited at 1000°c. 


4. The residue was fused with potassium bisulphate, dissolved in water, 
and a second tannin precipitation carried out. The precipitate was 
washed, dried, and ignited at 1000°c. 

5. The residue was fused with potassium bisulphate, the pH adjusted 
to approximately 4 by addition of sodium hydroxide solution. The iron 
was precipitated as hydroxide, leaving the aluminium in solution. This 
ferric hydroxide precipitate served to ‘ scavenge ’ radioactive impurities, 
especially uranium decay products, from solution. 


6. The pH of the solution was adjusted to approximately 4, a tannin 
precipitation carried out, the precipitate dried and ignited at 1000°c. 

7. The residue was fused with postassium bisulphate, and dissolved 
in water. A trace of uranium was still present in the solution. To the 
solution was added methyl orange indicator, followed by ammonium 
hydroxide solution, drop by drop, until the indicator turned yellow. 
The precipitate was centrifuged off, digested with ammonium carbonate 
solution (to remove uranium) and washed with water. 


8. At this stage the only activities detected with a gamma-ray spectro- 
meter were 7°Al, and, rather surprisingly "Na. The precipitate was 
dissolved in dilute nitric acid, sodium chloride ‘ hold-back carrier ’ added, 
and aluminium hydroxide precipitated by addition of 10°%, hexamine 
solution. 

9. The precipitate was washed with sodium chloride solution, taken up 
in nitric acid and the hexamine precipitation in presence of sodium 
‘hold-back carrier ’ repeated. 

10. The precipitate was taken up in a small volume of sodium hydroxide 
solution, a trace of solid residue was centrifuged off, dilute sodium chloride 
solution added, and aluminium hydroxide precipitated by addition of 
ammonium chloride. 

11. The precipitate was washed and dissolved in the smallest volume of 
dilute nitric acid. Examination of the solution with a gamma-ray 
spectrometer showed that the 1-28 mev peak of **Na was now absent. 
Ammonia was then added to the solution to precipitate aluminium 
hydroxide. 

When held at red heat on a platinum foil for several minutes the 
aluminium hydroxide decomposed to aluminium oxide, which formed a 
source suitable for use in the cloud chamber. 
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ABSTRACT 


The magnetizations of over 400 basalt specimens from the Deccan Traps 
have been found to be aligned in general along a NNW-SSE to NW-SE axis 
dipping towards the south. The results are consistent with the supposition 
that India has drifted northwards through over 50° of latitude and rotated 
25° anticlockwise within the last 70 million years. Im all three localities, 
the poles corresponding to the upper rocks are nearer the present geographical 
poles than those from the lower flows. The discrepancy may represent a 
northward drift of India during the formation of the Deccan Traps, but there 
are other possible explanations and further data are necessary to elucidate the 
problem. 


§ 1. LyTRODUCTION 


AN earlier survey of the Deccan Traps (Clegg et al. 1956) was undertaken 
by the Physics Department of Imperial College in the winter of 1954-55 
to initiate a project of palaeomagnetic research in India. Some 450 
basalt specimens from two areas, Khandala and Linga (fig. 1, (1) and (2)), 
were examined at that time and shown to possess on the whole consistent 
remanent magnetizations whose direction is almost perpendicular to that 
of the present earth’s field. The mean declination and dip for the two 
localities were N 155° E and 53° downwards respectively. This result ig 
in good agreement with that found independently by Irving (1956), and is 
most plausibly interpreted on the assumption that at the time when the 
rocks were laid down, during the late Cretaceous or early Eocene periods, 
India occupied a position considerably to the south of the equator, and 
that it has since drifted northwards and rotated in an anticlockwise 
direction. The ancient positions of the North and South poles corre- 
sponding to this direction of magnetization differed strikingly from those 
obtained by other workers for British, North American and Icelandic 
rocks of the same period, and this suggested that at least part of the 
postulated movement of India may have been in relation to other con- 
tinents. 

The present investigation has been carried out with the aim of extending 
the previous work to other areas in the Deccan Traps. This is of con- 
siderable importance, since the formation covers a large area in Western 
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India (about 200 000 square miles), and attains a thickness of several 
thousand feet in some places, and it is highly probable that rocks from 
different areas may differ considerably in age. Moreover, by extending 
the survey, it was possible to compare magnetic results obtained from 
rocks of various petrological types, for in a number of localities in the 
Deccan Traps acid volcanic rocks occur in addition to the much more 
prevalent plateau basalts. One of the best known of such localities is 
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Palaeomagnetic Surveys (Deccan Traps). 


(1) Khandala area. (5) Pavagadh area. 
(2) Linga area. (6) Kambatki area. 
(3) Nipani area. (7) Igatpuri Ghats, 
(4) Amba area, 
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Mt. Pavagadh (fig. | (5)), a prominent hill northeast of the city of Baroda 
in the north of Bombay Province where basic as well as acid lavas were 
collected. Inthe south of Bombay Province, 400 miles from Mt. Pavagadh, 
two other localities were sampled in detail : the Nipani area (3), near the 
small town of that name on the Kolhapur—Belgaum highway, and the 
Amba area (4), constituting a hill of the Western Ghats about 50 miles 
northwest of Nipani on the Kolhapur—Ratnagiri highway. Apart from 
the three major surveys, two small test collections involving a small 
number of samples were made from two hills, the Kambatki (6) and 
Igatpuri (7) Ghats, situated between the two southern areas and Mt. 
Pavagadh. <A preliminary and less complete account of measurements on 
rocks from Nipani, Amba, and Mt. Pavagadh has been given previously 


by Clegg et al. (1957). 


§ 2, COLLECTION AND MEASUREMENT 


Between October 1955 and November 1956, we undertook a number of 
field trips to the localities concerned and collected altogether 157 blocks of 
rock. At the collecting sites, these were marked with reference to the 
horizontal plane and to magnetic north, and the majority were taken to 
Bombay, where a palaeomagnetic laboratory suitable for the investigation 
of weakly magnetized geological specimens had been set up at the Tata 
Institute. Some of the rocks were shipped to Imperial College, London, 
for study. In each laboratory cylinders were cut from the blocks and the 
direction and intensity of their remanent magnetization were measured 
by means of a simplified adaptation of the astatic magnetometer described 
by Blackett (1952). In this, the standard procedure used by Clegg et al. 
(1954) was followed. In the present paper, the word ° specimen ’ will be 
used to denote each cylinder cut from a block, and the word ‘ sample ’ 
to denote a block. Altogether 480 specimens were measured, or an average 
of roughly three specimens to a sample. The results for the three areas 
are described separately below. 


§ 3. THE Nripant AREA 


The Deccan Trap flows in this region are exposed mainly on a number of 
low hills and in quarries. The present collection was made from two series 
of outcrops of basaltic lava, one occurring on the prominent escarpment 
lying three miles south of Nipani, and the other in two quarries on the 
plain ten miles north of Nipani. Samples were taken on the escarpment 
from five sites representing at least three distinct flows and a vertical 
column of nearly 200 feet. The single flow exposed in the quarries was 
650 feet below the highest flow sampled on the escarpment, but since the 
two localities were about 13 miles apart, a regional northward tilt of the 
flows by about 3° would suffice to place the quarry flow above the escarp- 
ment flows on the stratigraphic column. Despite the remarkable hori- 
zontality of the flows characteristic of most of the Deccan Traps (Wadia 
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1953), such a possibility cannot be completely ruled out. However, for the 


purpose of this paper, the quarries will be referred to as ‘lower Nipani’, 
and the escarpment as ‘ upper Nipani ’. 

The results are shown in table 1 and the directions of magnetization of 
individual specimens are plotted in fig. 2. Specimens cut from one 
sample in each of the lower and upper sections showed scattered directions 
of magnetization, and these have been omitted. It is evident from fig. 2 


Fig. 2 
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that the magnetic directions are grouped into two main classes, the 
ower Nipani specimens generally having SSE declinations with rather 
steep downward dips, while most of the upper Nipani specimens have 
NNW declinations with shallower upward dips. 


$4. THe Ampa AREA 


A vertical sequence of at least 1600 feet of basaltic lava flows is partly 
exposed on the steep side of the Amba Ghat, which lies on the western 
edge of the Deccan plateau just above the coastal plain. In the topo- 
oraphy of the Western Ghats this hill occupies a place corresponding to the 


Amba Area. 
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Bhor Ghat at Khandala, 130 miles to the north, mentioned previously by 
Clegg et al. (1956). Hight sites representing five or more distinct flows were 
sampled and again the polarizations fell into two major groups (fig. 3 and 
table 1), with most of the directions aligned approximately along a NNW- 
SSE axis dipping downwards towards the south, similar to the preferential 
direction at Nipani. A sub-division into ‘upper’ and ‘lower’ Amba 
areas was made to distinguish the three highest sites, which were mag- 
netized towards the NNW with an upward dip, from the five lower sites 
where the great majority of the samples had polarizations directed 
towards the SSE, but with rather steeper downward dips. As found in 
the Nipani area, a few of the lower Amba specimens yielded consistently 
odd results, with generally SSE declinations as expected, but also steep 
upward, instead of downward dips. Since 17 specimens from two samples 
and two different flows were here involved, values for their mean direction 
have been listed separately in table 1. 


§5. THe PavagapH AREA 


Mount Pavagadh is separated from the main body of the Deccan Traps 
by a distance of about 30 miles. It rises abruptly from the plains of 
Gajerat to a height of 2811 feet above sea level. The main mass of the 
hill, up to a height of about 2400 feet, consists of numerous nearly hori- 
zontally bedded flows of porphorytic as well as non-porphorytic basalt 
and olivine dolerite, together with some basic tuffs. Above this level the 
extensive Mauliya plateau is covered with green rhyolitic tuffs, and is 
capped with a large mass of rhyolite which rises steeply from the plateau 
to the summit. These acid rocks were first identified by Fermor (1906) 
who suggested that they might be interbedded with basic lavas, but 
Mathur and Dubey (1934) have since shown them to be younger than the 
basic lavas, and it is now believed that both the rhyolitic tuffs and the 
overlying rhyolites belong to an acid cycle of igneous activity which may 
have taken place considerably later than that which produced the basic 
rocks below. In the present account (see table 1) the lower basic rocks 
will be referred to as the ‘ Lower Pavagadh ’ and the upper acid tuffs and 
rhyolites will be termed ‘ Upper Pavagadh ’. 

Hight flows of largely basic composition were investigated in the lower 
Pavagadh area. The magnetic directions (fig. 4) were predominantly 
N to NNW declinations with shallow upward dips. Fifteen specimens 
from eight samples of acid tuffs of the upper Pavagadh area were also 
examined, Their magnetizations were found to be directed nearly 
towards the north, as in the basic rocks, and the dips were again shallow, 
but slightly downwards. The data for the tuffs are plotted as open 
circles in fig. 5. The mean direction of magnetization is 5° west of north’ 
with shallow downward dip. 

Directions for 69 specimens cut from 14 samples of rhyolite are also plotted 
on the same diagram. The results show too much scatter to permit the 
computation of a significant mean direction of magnetization. Different 
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specimens of rhyolite, even from samples collected at points not far apart, 
gave in some cases results which completely disagreed with one another, 
and they reacted quite differently to simple tests for their stability. The 
close spacing of the cluster of dots on the lower left side of fig. 5 would 
alone suggest that the 16 specimens concerned have stable magnetizations, 
particularly in view of the large angles which these make with the geo- 
magnetic field. To test this, a number of the specimens were subjected 
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to a.c. demagnetizing fields of 200 oersteds r.m.s.; this had no effect on 
their remanent polarizations. At the other extreme, when some of the 
specimens from samples showing much scatter were placed for periods 
varying from a few hours to a few weeks with their magnetization axes at 
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right angles to the earth’s field, and re-measured, this axis was usually 
found to have been pulled towards the magnetic meridian, in some cases 
through a complete right angle. 


§ 6. TEst CoLLECTIONS 


Four samples each from two flows of the Kambatki Ghat south of 
Poona, and from four flows of the Igatpuri Ghat northeast of Bombay, 
were examined. The lower Kambatki flow gave a mean declination of 
176° and a mean dip of +60°, whilst the corresponding values for the three 
upper Igatpuri flows were 161° and +51°. These results are similar to 
those obtained for the Khandala, Linga, lower Nipani and lower Amba 
areas, but not enough samples have been examined to justify any further 
conclusions. The upper flow at the Kambatki Ghat and the lowest flow 
at the Igatpuri Ghat had respectively NE and NW declinations with steep 
downward dips, perhaps due to partial instability. 


§ 7. Discussion oF RESULTS 


The complete results for the Nipani, Amba, and Pavagadh areas are 
summarized in table 1. The general direction of magnetization is along 
a NW-SE to NNW-SSE axis dipping towards the south. This is in 
broad agreement with the earlier results of Irving (1956) and Clegg e¢ al. 
(1956), and supports the supposition that India has moved northwards 
across the Equator relative to the earth’s geographical axis, and rotated in 
an anticlockwise direction, since the rocks were laid down. In the 
Nipani and Amba areas the lower flows have reverse (southeasterly) 
magnetizations with rather steep downward dips, while the uppermost 
flows are normally magnetized and the dips are significantly less. In the 
Pavagadh area there is also a significant difference between the lower 
basic rocks and the acid tuffs above, both having normal polarizations 
and shallow dips, but whereas the dips of the lower flows are upwards, 
those of the overlying tuffs are downwards. No significant mean direction 
could be computed for the rhyolities overlying the tuffs. 


7.1. The Directions of the Magnetic Axes and Ancient Pole Positions 


Since the period of time over which the rocks from each individual site 
were extruded is likely to have been large in comparison with the period 
of secular variation, it is reasonable to suppose that the magnetic directions 
shown in table 1 correspond to those of a central magnetic dipole coaxial 
with the earth’s rotational axis. In column 4 of table 2 the values of 
the ancient latitude A of the various sites have been computed on this 
assumption, from the mean dip 0, using the formula 


tan A= tan 0. 
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The ‘ reversely’ and ‘normally’ magnetized series of flows have been 
grouped together for comparison. Column 5 of table 2 shows the ancient 
positions of the North geographical pole calculated from the values of A, 
and from the mean declinations. These pole positions are compared with 
those obtained previously by Clegg et al. (1956). The few anomalous 
results from the various sites (e.g. the 17 specimens from the lower Amba 
area showing upward dips) have been omitted from this table. 

The discrepancies between the values of A shown in column 4 arise 
partially because the various sites are separated by distances of several 
hundreds of miles, and a better comparison is obtained by referring all 
the results to a single locality. For this purpose Bombay (present latitude 
19° N) has been chosen, and columns 6 and 7 of table 2 show the ancient 
latitude and the ancient axial dipole field dip at Bombay, corresponding 
to the pole positions in column 5. 


7.2. Possible Implications of the Results 


From the summarized results set out in table 2 it is apparent that the 
rocks from the various localities visited fall into three main groups, viz.: 


(a) Those of the lower Nipani and lower Amba sections, for which the 
magnetic directions and ancient pole positions correspond reasonably well 
to those obtained previously for the Linga and Khandala areas by Clegg 
et al. (1956). 

(6) Those of the upper Nipani and upper Amba sections together with 
the basic rocks of lower Pavagadh, which are normally magnetized and 
have markedly shallower upward dips, and for which the corresponding 
ancient pole positions lie considerably closer to the present geographical 
pole. 

(c) The upper rhyolitic tuffs of upper Pavagadh, which are again 
normally magnetized, but which have a shallow downward dip corre- 
sponding to a low northerly latitude. 


The mean directions of magnetization for these three groups, obtained 
from table 2 (directions are referred to Bombay) are 


Mean declination Mean dip 
(a) N 156° E +56° down 
(bd) N 341° E — 24° up 
(c) N 355° EH +10° down 


One obvious possible explanation of these facts is that the three groups 
of rocks were laid down during different geological epochs, between each 
of which India moved northward relative to the earth’s geographical 
axis. It may also be that the earth’s magnetic field underwent a reversal 
between the first and second of the epochs. The geological evidence 
provides some support for this view. Thus in the Amba Area there is a 
section 200-300 feet thick between the uppermost ‘ reversed ’ and lower 
‘normal’ flows in which no exposures were found, so that the upper and 
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lower sections probably differ considerably from one another in age. 
Again in the Nipani area a horizontal gap of 13 miles, containing four 
outcrops, separates the two localities, so that a strong possibility exists 
that the missing section represents a relatively long time interval, although 
it is not definitely established that the ‘ upper’ Nipani rocks are indeed — 
younger than those of the ‘ lower’ Nipani area. The unusual sequence on 
Mt. Pavagadh is of particular interest. The lower basic rocks may well 
belong to the uppermost flows of the Deccan Traps, and it is probable 
that a considerable time elapsed between the last basalt outpouring and 
the beginning of the acid cycle responsible for the overlying rhyolitic 


Fig. 6 


I. Pavagadh (Upper) V. Linga 

II. Pavagadh (Lower) VI. Nipani (Lower) 
IIT. Amba (Upper) VII. Amba (Lower) 
IV. Nipani (Upper) VIII. Khandala 


tufts. Mathur and Dubey (1934) quote results of radioactive age deter- 
minations of basalt and rhyolite specimens from Mt. Pavagadh, made by 
the helium method. Although helium dating seldom yields precise 
results they estimate from their data that the age of the basalt is early 
Kocene, in conformity with the remainder of the Deccan Traps, while 
that of the rhyolite is Miocene. If we assume the position of ane at the 
end of the Cretaceous period to be that originally postulated by Clegg 
et al. (1956) and by Irving (1956), and the drift to have proceeded at a more 
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or less constant rate up to the present time, the area under investigation 
should have crossed the equator during the early Miocene. The observed 
mean dip of the tuffs would tend to support this supposition. 

On the other hand, it is by no means certain, on the basis of the evidence 
so far available, that this variation of magnetic dip with depth can be 
taken as positive evidence of a continuous northward movement of India. 
One major difficulty in the interpretation of the results arises from the 
uncertainty in the relative dating of the basic rocks from the various sites. 
For the total depth of the Deccan Traps, particularly towards their eastern 
edge, is known to exceed 10 000 feet in places (Wadia 1953), and as it was 
only possible to sample a small fraction of this range at any of the sites 
visited, it seems improbable on primae faciae considerations that the rocks 

of either of the two categories (a) and (b) above are from corresponding 
parts of the stratigraphical column. 

There are obviously other possible physical causes which may have given 
rise to the difference in inclination between the categories (a), (b) and (c). 
It may, for instance, be due to some physical phenomenon associated with 
different depths of burial, such as the magnetostriction effect suggested by 
Graham (1957). Alternatively it may result from partial magnetic insta- 
bility, and the consequent pulling of the magnetic directions towards that 
of the present earth’s field. For rocks which are magnetized normally 
this latter effect would tend to reduce the magnetic dip, while for those 
magnetized reversely it would tend to increase the dip. In this connec- 
tion it is worth noting that the rocks of category (a) which have the 
steepest dips, are reversely magnetized, while those of categories (6) and 
(c), with shallower dips, are all normal. 

Work is continuing in an attempt to find exposures of intermediate 
sections in the areas sampled so far, and the survey is also being extended 
to include a detailed study of vertical sequences in other parts of the 
Deccan Traps. 
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ABSTRACT 


Three distinct lines of experimental work now point towards a new con- 
ception of the electronic structure of the transition metals. Recent x-ray 
work has shown that in Fe and Cr the valence electrons are not as well 
localized near the atoms as had been supposed, while in Co, Ni and Cu tho 
compact 3d-shell is nearly complete. Saturation magnetic moments of Ni 
rich alloys show that Ti and V act as though they had respectively four 
and five electrons to contribute to the conduction band of Ni. Neutron 
diffraction evidence shows that pure Ti and V have no magnetic moments, 
ordered or disordered, and that in Ni-Cr alloys the Cr atom carries a moment 
not much different from Ni, though each added atom removes nearly 5 uB 
from the saturation moment. We therefore suggest that elements of atomic 
nuinber less than Fe have few 3d electrons and many conduction electrons, 
while for those above iron this is reversed. It appears as though this 
conclusion can be maintained almost independently of crystal structure or 
composition of alloy. On this basis the variation of saturation magnetic 
moment of all the binary ferromagnetic alloys can be explained. 


$1. INTRODUCTION 


THE electronic structure of the first long period elements has been 
discussed for many years (see e.g. Hume-Rothery and Coles 1954). 
In particular the ferromagnetic elements have been of interest, and 
there has been much debate about the most suitable model for their 
description. Recently the technique of neutron diffraction has allowed 
us to learn a great deal more about the spatial distribution of magnetic 
moment in both pure metals and alloys. In addition, a recent x-ray 
study by Weiss (1957) has shown that the spatial distribution of electron 
density is very different in the metals Cr and Fe from what it is in the 
free atoms of these elements. In this paper we outline a scheme of 
metallic valences which is directly suggested by the new experimental 
results and which appears to correlate a wide range of observations 
with a minimum of assumptions. In § 2 we outline the experimental 
results, and in §3 state the main assumptions and implications of 
the model proposed. § 4 is devoted to a discussion of the pure metals 
from vanadium to zinc. §5 handles the binary magnetic alloys. $6 
includes some speculations on the interrelation of electronic structure 
.and crystal structure. 
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§ 2. EXPERIMENTAL RESULTS 
2.1. X-ray Electron Counts 


Weiss (1957) reports the results of careful determination of the 


crystal form factors for the pure metals Cr, Fe, Co, Ni, Cu. He separates 
the electron density into three parts ; the compact argon core containing 
18 electrons, the electrons whose radial density distribution resembles 
atomic 3d functions, and those roughly uniformly distributed over 
the crystal. The contribution of the argon core can be accurately 
calculated, and Weiss discusses the residual scattering power, which 1s 
to be interpreted as arising only from the d-like electron density. He 
measures the decrease in integrated intensity of reflection from certain 
planes as the scattering angle is increased. This decrease is interpreted 
as due to the primary extinction, and is consistent with the d-like radial 
distribution. The absolute intensity extrapolated to zero scattering 
angle leads to the following numbers of electrons in d-like orbitals : 
Cr, 0-2; Fe, 2:3; Co, 8-4; Ni, 9-5; Cu, 10. 

The accuracy is quoted to be -+0-3 electron. The remaining valence 
electrons are concluded to be in conduction states, with a roughly 
uniform density distribution. The number of conduction electrons are 
therefore about: Cr, 5:8; Fe, 5-7; Co, 0-6; Ni, 0-5; Cu,1. While 
the experiment says little about the actual distribution of these electrons, 
they certainly do not show a marked density peak in the same position 
as that of the atomic 3d functions. For convenience we shall denote 
these states as 4c states, indicating that they are conduction band states 
built up mainly from 4s and 4p atomic states (but in fact containing 
some admixture of all higher states). 


2.2. Neutron Diffraction 

It is found (Shull and Wilkinson 1955) that the neutron diffraction 
from ferromagnetic crystals is consistent with a distribution of electron 
spin corresponding to 3d functions around each atom. The total electron 
spin so measured is consistent with the observed saturation magnetic 
moment, so that we may conclude that the magnetic electrons are 
3d-like in electron density ; for Fe, Co and Ni there are about 2, 13, 4 
respectively. 

Above the Curie point (Wilkinson and Shull 1956, Lowde 1956) spin 
disorder scattering is found. Although the mean magnetic moment is 
zero, it appears that every atom carries the same spin as in the magnetized 
state, but now randomly oriented except for a certain amount of short 
range order (Elliott and Marshall 1957). To describe such a system, in 
which the spin on each atom is well-defined, the Van Vleck (1953) model 
seems appropriate. 

In ferromagnetic alloys, neutron observations demonstrate that in 
general terms each atom tends to carry approximately the moment 
appropriate to it in the pure metal (Shull and Wilkinson 1955). Each 
alloy system really merits detailed discussion (§ 5). 


~ 
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2.3. Magnetic Measurements 


The course of the saturation ferromagnetic moment. in many alloy 
systems has been studied. The results are summarized by (e.g. Bozorth 
1951), but in many instances it is necessary to consult the original papers 
for detail. Many of the alloys follow the Slater—Pauling curve (fig. 1), 
but many do not. The strongest deviations from this curve, in the 
Ni-Cr, Ni-V, and Co-Cr series of alloys, will receive a ready explanation 
on the proposed scheme. 


Fig. 1 
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Slater—Pauling curve of mean saturation magnetic moment per atom in binary 
alloy systems plotted against number of electrons per atom (after 
Bozorth). 


§ 3. THE PROPOSED SCHEME OF STRUCTURES 


The proposed scheme is based on the separation of the electronic 
states into 3d-like states and conduction states, the number of electrons 
in the 3d-like states for the pure metals being just as found by Weiss. 
It is by no means easy to determine theoretically what a full self- 
consistent field solution for any particular metal would be like, and 
information on free atoms is of little direct use. Our picture is that 
as the nuclear charge is increased and electrons are added to the argon 
core from potassium right through to chromium, there are practically no 
3d-like electrons; they are all ‘conduction electrons’ with spatial 
density distribution quite unlike that of the atomic 3d-states. At 
chromium the d-levels begin to fill. At cobalt there is a sudden change 
to a configuration with a high number of d-like electrons and few conduc- 
tion electrons. This type of structure of course persists through copper, 
zinc and gallium with their intact (3d)}° filled shells. We make no 
attempt at the moment to give this picture theoretical backing. Note 
that the main determining feature of the metal which fixes the number 
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of d electrons is the nuclear charge of the atom; the crystal structure 
we would choose to regard as a consequence of the electronic structure, 
playing only a minor role in determining just what that electronic 
structure is to be. 

Now let us take a brief look at the expectations for alloys. Here we use 
the model suggested by Friedel (1954), which reconciles the localization of 
electron charge near positively charged impurities with the original ideas 
of Hume-Rothery ef al. (1934) and Jones (1934) of adding electrons 
to a nearly rigid conduction electron band structure. This is to be 
considered along with the existence of an equilibrium between the 
conduction and d state occupations. In alloying elements of similar 
electronic structure, i.e. high d-like elements amongst themselves or 
low d-like elements amongst themselves we expect only smooth variations 
of macroscopic quantities such as the mean magnetic moment—although, 
as we shall see in § 5, we do not expect this variation to be of the same 
type in the two cases. But now let us examine the case of a high d 
element dissolved in a low d one, e.g. Niin Fe. The conduction electrons 
cannot be localized very effectively in most metals, so that a nickel atom 
substituted into an iron lattice must find considerable 4c electron density 
overlapping onto it. The only reasonable recourse open to restore 
approximate electrical neutrality in the cell occupied by the nickel atom 
is for it to delocalize some of its d-like electrons, by promoting them 
into the conduction band of the alloy. Thus the high d elements may 
lose some d electrons when dissolved in low d elements. 


———— ee 


The remaining case to consider is that of a low d element, say V, — 


dissolved in a high d element, say Ni. The five 4c electrons of the 
vanadium may be regarded as added to the nickel conduction band. They 
therefore disturb the equilibrium between conduction and d-like levels, 
with consequent changes in magnetic properties. Friedel’s description 
in terms of the V atom removing localized states from the conduction 
band is preferable to the more naive concept of adding conduction 
electrons to the conduction band, and then allowing them to overspill 
into the high-density states provided by the d band. Whichever language 
we use, the essence of our suggestion is that V is five times as effective 
as Cu in increasing the occupation of d states in Ni. 
Let us then summarize our hypotheses. 


(1) The valence electrons may be divided into two classes—conduction 
electrons distributed over the metal, and d electrons concentrated near 
each atom. 


(2) The magnetic properties are due essentially to the spin moment 
of d electrons only. 


(3) The interactions between magnetic moments of near-neighbour 
atoms in a crystal are unknown even as to sign; minor component 
moments may be parallel or antiparallel to a ferromagnetic solvent 
moment. 
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(4) For some metals, the number of d electrons is an integer, while 
for others two different configurations may be in dynamic equilibrium 
with the conduction band (Mott and Stevens 1957) 


(5) In alloys of whose components have roughly the same conduction 
electron density, each atom retains approximately the d electron con- 
figuration characteristic of the pure metal. 


(6) In alloys where the conduction electron densities of the pure 
metals are very different, the number of d electrons on the minor 
component may change to decrease this difference. Mn and Fe in 
high-d solvents become high-d atoms ; Ni and Co in low d solvents lose 
one or two d electrons (starred columns in table 1). 


(7) In alloys whose major component has equilibrium between con- 
duction and d states, each addition of a conduction electron by alloying 
causes the occupation of one d state. 


§ 4. THe Pure Merats 


Cu. The electronic structure of metallic copper can scarcely be 
questioned. We concur with the usual assignment (3d)! (4s)!. 

Ni. The structure of nickel lies between (3d)!° and (3d)9 (4s)!. The 
precise description. adopted depends on the view we take of spatial 
correlations of d-like and conduction electrons. If all correlation is 
neglected we would arrive at the simplest form of ‘ collective electron ’ 
treatment in terms of ordinary d and s bands. However, we prefer to 
believe that the number of (3d)° ions is very low—i.e. that two 3d holes 
are seldom found on the same atom. We therefore think of the d 
electrons in terms of the Van Vleck model, in which the proportions of 
(3d)!° and (3d)° configurations are in dynamic equilibrium, and in which 
(3d)* configurations do not occurt. Above the Curie point the spin 
orientation of the (3d) configurations shows no long range order, and so 
gives rise to diffuse spin disorder neutron scattering. For pure nickel 
the proportion of (3d)° is about 60°, and (3d)!° about 40% so that 
the mean spin moment per atom is 0-6 wp, as observed. It may be that 
in fact the proportion of (3d)® is somewhat lower than 60°, for there will 
be some contribution to the observed moment both from imperfectly 
quenched orbital magnetic moment, and from partial spin alignment of 
the conduction electrons. 

We regard the 4s electrons (0-6 per atom) as uncorrelated with the 
d electron holes, and so use a simple band of conduction states. 

Co. The discussion for cobalt is essentially similar; the two con- 
figurations involved now are (3d)° and (3d)* (4s) carrying spin moments 


+The actual proportion of (3d)* states in pure nickel is not known. It 
simplifies our discussion of saturation magnetization to suppose it 1s low. It is 
hoped to discuss the value of the Curie temperature using this supposition in a 
later paper. The general picture of the alloying behaviour of saturation 
moments is not, however, sensitive to the supposition, 
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of one and two Bohr magnetons respectively. Neglecting orbital and 
conduction electron contributions to the saturation magnetization the 
value fi=1-7 yy for Co indicates 70%, (3d)* and 30%, *(s0)* 

Fe. The iron atom we believe may exist in two distinct states; a 
high-d state analogous to Ni and Co (occurring in some alloys), or in a 
low-d state. Body-centred cubic iron we think is the latter, being in the 
pure configuration (3d)? (4c)® There is no equilibrium between this 
and other states, and the d electrons do not contribute to the Fermi 
surface (Mott and Stevens 1957). The two aligned electron spins contri- 
bute an effective magnetic moment of 2,. Measurements of the 
gyromagnetic ratio show that the orbital momentum contribution is 
about O-lj,. The remainder of the observed saturation moment, 
2-22 u,/atom, may be contributed by conduction electron polarization. 
Griffith (1956) has suggested that the d orbitals concerned are the Kg 
orbitals which are split off from the rest by the crystalline field. This 
assignment will allow us to assume that the orbital magnetic moment 
is well quenched as is observed. 

Mn. Manganese is probably similar to iron. The pure metal may 
contain both (3d)! (4c)*® and (3d)? (4c)? configurations in dynamic equili- 
brium, and in some alloys Mn may take up configurations (3d)? or (3d)° 4s. 

Cr. Chromium appears to be a mixture of (3d)! (4c)? and (4c)® con- 
figurations. The antiferromagnetic structure (Shull and Wilkinson 1953) 
is attributed to the spin of the former configuration. 

V. Vanadium is probably pure (4c)®> configurations, for Shull finds 
no magnetic scattering (Shull and Wilkinson 1953). 

Ti. Titanium is probably pure (4c)* on the same reasoning. 


§5. THe FERROMAGNETIC BINARY ALLOY SYSTEMS 


In this section we apply the ideas of the preceding sections to the 
saturation magnetization curves of binary alloys, and to the moment 
carried on the atoms of each constituent ; results of experiments are 
available for saturation moments of most systems, but only a few neutron 
diffraction results are available to apportion the total moment between 
constituents. For each dilute alloy we shall try to assign an electronic 
configuration to both matrix and solute. It is satisfactory that treating 
each alloy on its own merits in this empirical way leads to a consistent 
and plausible scheme of structures, summarized in table 1. 

Ideally, we want the results for very dilute alloys, so that we may 
neglect interaction between solute atoms. These results are of course the 
most difficult to get. Only in exceptional cases, where a linear extra- 
polation of results at moderate concentration back to zero concentration is 
acceptable, can the variation of mean moment be given to better than 
++0-2 4g per solute atom. In systems with limited primary solubility 
ranges (notably Fe Mn, Co Fe, Fe Ni, Co Al) the uncertainty is considerably 
greater. 

It is also true that it is not possible to assign an accurate magnetic 
moment to an atom of a given electronic structure ; to a first approximation 
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we may suppose that the 3d electrons always produce a state of 
maximum possible spin, and that this is the whole observed magnetic 
moment. But in fact terms of the order of 0-1 or 0-24, per atom may 
arise from incompletely quenched orbital motion and from spin polariza- 
tion of the conduction band electrons. We cannot therefore expect 
perfect agreement with experiment; discrepancies of up to 0-2 or 0:3 pr, 
between assignment and observation are not immediately relevant. 


Table 1 
Ti V Cr Mn Mn*| Fe Fe* 
P fees nee 8 | ss CeO. 6 sh et 0. 
d 0 0 i () Wi es Gl i 2 0 as 2 
je 0 0 1a) fy eS Gl 3 2 oe 2 
4 low d elements 
SS Se es cm 
Co Co* Ni Ni* Cu | Zr 
Cc i () a a 1 <> 0 Oo as Hy i) 2 
d fs) <es WY) () ss 7/ oF 10 —— 10 10 
je Ae ae | Gens 3 1 (0) of? 0 0 
high d elements ——HW¥—___—____+> 


The number of conduction electrons c, number of d electrons d, and atomic 
magnetic moment 1 for atoms of the first transition period. The starred columns 
refer to the behaviour of the element in solutes of the opposite type con- 
figuration. Where no starred column is given the behaviour is the same in both 
high and low-d solvents. Equilibrium between structures is indicated by 
horizontal arrows. 


The Slater—Pauling curve of average saturation magnetic moment of the 
transition metal alloys provides an invaluable general background to our 
discussion. The curve is shown in fig. | which plots the known behaviour 
of many alloys against their effective electron configuration. The alloys 
of adjacent metals define the course of the main curve. Roughly speaking 
we see that alloys of low d elements with iron form branch A of the curve, 
while alloys of high d elements with other high d elements form branch B. 
The maximum of the curve, part C, is formed from alloys of iron with 
high delements. Alloys of low d elements other than iron with the high d 
elements, cobalt and nickel, give the curves D showing very strong 
deviations from the main curves. We shall now discuss these four types 
of alloys A, B, Cand Din turn. A summary of this discussion is given in 

‘table 2, 
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5.1. Type A Alloys 


These are alloys of low-d elements, Al, Si, Ti, V, Cr and Mn with Fe, 
and should behave according to hypothesis (5), §3. Because Fe is in a 
pure (3d)2(4c)® configuration in which the magnetic electrons do not 
contribute to the Fermi surface, replacing an iron atom by any elenent 
does not affect the magnetic properties of the other iron atoms. The 
mean magnetic moment will decrease by 2-2, for each magnetically 
inert atom added. Cr and Mn contribute their own 3d spin moments, 
either antiparallel or parallel to the iron 3d spins, giving correspondingly 
more or less rapid decreases of saturation moment. Let us consider the 
alloys one by one. 

Fe Al. The initial slope of the moment/concentration curve is observed 
to be —2:3 (Fallot 1936), in good agreement with the expected — 2-2. 
At 18%, Al superlattice lines appear and the moment falls more rapidly. 

Fe Si. The initial slope is — 2-2 (Fallot 1936). The moment falls more 
rapidly beyond 7° Si, where superlattice lines appear. 

Fe Sc, Fe Ti, have not been studied ; we expect a slope of —2:-2. 

Fe-V. Experimentally (Fallot 1936) the slope is observed to be — 2:2, 
as expected ; this holds up to 30°, V at which there is a phase change. 

Fe-Cr. Cr is a mixture of (4c)® and (4c)9(3d)! configurations. When 
Cr is added to iron, if the proportion of (4c)°(3d)! atoms is 2, the mean 
moment will vary by —2-2 —wz if the Cr spin is antiparallel to the iron, 
or —2-2 +z if the Cr spin is parallel to the iron spin. The experiments 
of Shull and Wilkinson (1955) show that the moment on Cr differs from 
that on Fe by -:3 wp in these alloys. We accept the solution that the 
Cr moment is about —1 pp, i.e., antiparallel to the iron moment. (The 
second solution fig,;=5 would require five d-like electrons, which does not 
fit our general scheme. In addition it implies a rapid variation of the 
moment on iron atoms, which is not seen in other systems.) The initial 
rate of decrease of saturation moment estimated from Fallot’s (1936) 
curves is —2-9. This suggests the Cr atoms enter 70°, as (4c)°(3d)!, and 
30% as (4c)®. This would not be inconsistent with extrapolation of 
Shull’s results (see fig. 2 (@)). 

At higher concentration of Cr, the experiments show the mean slope 
to be about —2-2 (the value quoted by Stoner 1947). The change in 
slope may be due to an antiferromagnetic coupling between neighbouring 
chromium atoms. If this coupling is more important than that between 
chromium and iron atoms, we would suppose that by the time a concentra- 
tion of about 40°, Cr is reached the net moment on Cr atoms will be zero, 
so that the mean slope over a wide concentration range corresponds to 
dilution of the iron, as observed. Shull and Wilkinson, however, shows 
a small positive moment on Cr atoms at 45°. It would be interesting to 
have confirmation of this, and to know if, say at 60°/ Cr the moment is still 
larger, for we would expect it to remain essentially zero. At still larger 
concentrations of Cr the antiferromagnetic structure of pure Cr must 
develop, and dissolved iron will be unaligned, so the mean moment on 
iron atoms must fall in Cr-rich alloys, as observed. 
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Fe-Mn. Mn we expect to enter Fe as a mixture of (3d)4(4c)® and 
(3d)2(4c)5 configurations. If the relative proportions are (l—a) to a, 
the small concentration results would be 

Pyn= + (1+2) 
dji/de= — 2-2+(1+-2). 

Sadron’s (1932) experimental results are completely at variance with this. 
He quotes dj/dc——2-0 at zero concentration. However, the solubility 
of Mn is only 3-5 atm’, in «-iron, so that it is difficult to obtain accurate 
enough results on solutions of accurately determined concentration. 
A new measurement for this system would be very welcome ; if Sadron’s 
result is confirmed it will be difficult to fit smoothly into our scheme, 
though it is possible that Mn enters iron as (4c)’. This would be surprising 
in view of the very well established result that magnetically active d 
electrons are present round Cr atoms dissolved in iron. 


5.2. Type B Alloys 


These are alloys of Ni and Co with other high d elements and hypothesis 
(5) is to be combined with hypothesis (7). The mean moment of Ni and 
Co atoms in these alloys is controlled by the equilibrium between the 
number of conduction electrons and the number of d holes. As indicated 
in § 2 the perturbation of the conduction band on adding alloying elements 
is equivalent to adding electrons to a rigid conduction band. But if these 
conduction states are in equilibrium with a set of d states with a high 
density of states, the net effect is the removal of a corresponding number 
of d-holes. This argument is independent of the use of conventional band 
or a Van Vleck description of the d-states, though they differ in their 
description of the distribution of spin moment in the alloys. 

Ni-Cu. Cu enters as (3d)! 4c, that is, carrying 0-4 more conduction 
electron than pure Ni giving a corresponding decrease in total Ni moment 
of 0-4 4,. The inert (3d)! replaces the mean Ni configuration of moment 
0-6 wp, and the total change is therefore — | w, for each added copper atom, 
just as observed (Marian 1937). On this model j,,=0, 
fyi = (0-6—c)/(1—c); on the simplest band theory jig, =fiyj=(0-6—c). 

Ni-Zn. Zn will be twice as effective as Cu in reducing the moment and 
so dju/dc should be —2. Experimentally it is observed to be —2-1 (Marian 
1937). 

Co-Cu. In this case the Cu replaces on average the configuration 
(3d)*3(4s)""?, 1.e., a moment of 1-7, is lost on average, and a further 
0-3 wp is lost from the localization of a further 0-3 conduction electrons 


around Cu. Hence da/dec=—2-0. There are no experimental results to 
compare to the prediction. 
Co-Zn. Here we expect dji/dc=—3-0 but no experimental results seem 


to be available. 
Ni-Co. Both pure metals are ferromagnets with approximately the 
same number of conduction electrons per atom. We therefore expect an 


~~ ants 
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almost linear variation of the mean moment with Co concentration from 
the initial value of 0-604 to the final value 1-714. This is experimentally 
observed; there is a slight change of slope at the f.c.c. > h.c.p. phase 
change. (Weiss and Forrer 1929). 


5.3. Type C Alloys 

These are alloys of Co and Ni with low 3d elements and show the 
greatest deviation from the Slater—Pauling curve. It has been suggested 
that the very rapid fall in the mean moment with increasing concentration 
of the low 3d element is due to large moments coupled antiparallel to the 
cobalt or nickel moments. But Shull and Wilkinson’s results on Co—Cr 
show that this is certainly not the case in this alloy but that the fall is due 
to the decrease in the mean Co moment on alloying. It therefore seems 
likely that in none of these alloys is there a large antiparallel moment on the 
low 3d element. A natural explanation is that the large number of con- 
duction electrons which these elements have enter the conduction band 
of the metal, disturbing the equilibrium between the various 3d configura- 
tions in nickel or cobalt and so suppressing the mean moment rapidly, 
according to hypothesis (7). 

Ni-Al. We expect dj/dc to be about —3-0. Experimentally (Marian 
1937) it is observed to be —3-0-+0-2. 

Ni-Si. Here we expect dj/dc to be —4:0. Marian (1937) has found 
experimentally a value —3-4. 

Ni-Ti. We suppose that titanium enters in the configuration (4c)* so 
we expect about dji/dc=—4. Experimentally the value is —3-8 (Marian 
1937). 

Ni-V. Vanadium has the configuration (4c)?so we expect approximately 
dj/dec=—5. Experimentally the value is —5-2 (Marian 1937). 

Ni-Cr. Cr enters Ni as a mixture of 3d1(4c)° and (4c)®. Experimen- 
tally (Marian 1937) dj/dc=—4-5, which is fitted if the weights of these 
configurations are respectively 75°%, to 25°%, with the spin on Cr atoms 
parallel to the Ni spins, giving fig, =0°75, fix; = 0-6-4650, dj/de=—4-5. 

Co—Al. Here we expect dj/de to be about —4-0. Experimentally 
Farcas (1937) has found a value of —2-4. This is a difficult alloy to work 
with as only 1°% of Al is soluble in Co at 300°c so this experimental result 
may be open to some doubt. 

Co-Si. We expect a value for dj/de about —5-0. Experimentally 
Sadron (1932) observed about —5-3. 

Co-Ti. We expect a value of about —5-0 but this alloy has not been 
examined experimentally yet. 

Co-V. We expect a value of about —6-0 but no experimental results 
are available. 

Co-Cr. Cr enters Co as a mixture of 3d4(4c)> and (4c)*; to fit the 
experimental results (Farcas 1937) the probability of (3d)! (4c)° is 0-3 and 
of (4c)® is 0:7 to give, for small concentrations, ji,,=0'3, Pg9=1-7-5¢ 
and dji/dc=— 6-4. 
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5.4. Type D Alloys 


These are the alloys which form the maximum of the Slater—Pauling 
curve and the explanation of their properties rests on hypothesis (6). 
It is convenient to treat not only Fe—Co and Fe—-Ni here, but to take 
all Fe and Mn alloys with high-d metals. There are two main ideas in 
this discussion. First, as described in § 3, high-d elements Ni and 
Co may be expected to lose one or two 3d electrons when they are 
dissolved in iron, in order to compensate the charge-overlap effect. 
The second assumption is that Fe and Mn are finely balanced between 
low-d and high-d so that when they dissolve in a high-d solvent the 
inverse of the previous effect causes them to gain not one or two d electrons, 
but to switch right over to a high-d configuration with less than one 
conduction electron. 

An unsatisfactory feature of the discussion of these alloys stems 
from our inability to predict the extent of quenching of orbital moment. 
A static cubic crystalline field quenches the orbital moment of (d)° 
and (d)® configurations, and the observed gyromagnetic ratios for Ni 
and Co confirm the smallness of orbital contribution for these configura- 
tions. But such static fields do not quench the orbital moment in (d)® 
and (d)? configurations, and we have no direct experimental information 
on the matter. It seems doubtful that we should really expect a large 
orbital contribution ; as a working hypothesis, we suppose that the 
spin moment is the major contribution for these configurations also. 

Fe-Co. Pure cobalt is a mixture of (3d)° and (3d)° configurations. 
The charge-overlap effect promotes some of these d electrons into the 
conduction band when Co is dissolved in iron; it may enter iron as 
either a mixture of (3d)® and (3d)’, or (38d)? and (3d). It is impossible 
to decide between these on the basis of the experiment, because we do 
not know how well the orbital momentum is quenched in the (d)* and 
(d)’ configurations. The spin moment of (3d)* alone would give 
Ho o=8 and da/dec—0-8. The experimental value (Weiss and Forrer 1929) 
of dujde is 1-2. If the orbital contribution of (d)’ is greater than 0-4 
there must be some admixture of (d)* configuration ; if it is less than 
0-4 there is some admixture of (d)°. 

Now consider cobalt-rich alloys. We suppose that small concentra- 
tions of iron dissolve in a high-d state, some mixture of (3d)? and (3d)8 
or possibly (3d)° and (3d). It is again impossible to assign weights 
to the mixture, for the orbital contributions are unknown. However, 
it is reasonable to say as this is now an alloy of type B, we expect Co Fe 
to form an extension of leg B of the Slater—Pauling curve. If the mixture 
of (d)’ and (d)* is considered, with weight « and 1—« respectively, and 
moments 3+-6, and 2+6, respectively, we get 


Roa tec Usi a) 
fepe=2+6,+2(1+6,—8,) 
djjdc=1+6,+2(6,—6,). 
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Here ¢ is the concentration of Fe atoms and the formulae are correct 
only to first order in c. Experimentally dji/de is observed to be 1-1 
so we can get agreement even with small values of 8,. (The observation 
is extremely uncertain in the hexagonal cobalt phase—the slope quoted 
is for f.c.c. phase (Weiss and Forrer 1929). 

Because we suppose that in iron rich alloys Fe behaves as a low-3d 
element but in cobalt rich alloys Fe behaves as a high-3d element it 
is clear that at some composition range a change over from one type to 
the other must occur. One possibility is that the change over occurs 
at the transition between the body-centred and face-centred cubic 
phases, i.e. at a concentration of about 75°, Co. However, we believe 
it to be more probable that the change over occurs at the maximum of 
the moment/concentration curve, i.e. at about 40°/, Co, as indicated in 
fig. 2 (c). (Shull and Wilkinson has found that in ordered Fe Co the Fe 
moment is 2-Syup. If at this concentration iron were still behaving 
as a low 3d element and only the Eg orbitals are available then the maxi- 
mum value for j.,,, would be about 2-2.) 

Fe—Ni. The discussion of this alloy is similar to that of Fe-Co. Pure 
nickel is a mixture of the configurations (3d)° and (3d)! Because of 
the change overlap effect in iron rich alloys it will have less 3d electrons 
and we suppose therefore that it goes in as a mixture of the configurations 
(3d)? and (3d)%. If the probability of the former configuration is x and 
of the latter (1—x) then 

fing =2+8 g-+2(1-+-8,—3,) 
[epee 
dji/de= —0-2+6,+a(1+6,—56,). 
The experimentally observed value of dji/dc is 0-6 so assuming 6,~0-2 we 
- get x~0-6/0-8+5,). 

For the nickel rich alloys we suppose that iron goes in as a mixture 
of (3d)? with probability y, and (3d)® with probability (l1—y). Then 
jigy=2-+8 9X y(1+82—85) 

fe — 0-6 c(0-6—y) 
djidc=2+8,+-y(07—9s). 
Shull and Wilkinson (1955) results on Ni-Fe show that jy, is roughly 
constant so we can set y=0-6. Experiments indicate that at zero 
concentration ji,,.223-0, and dji/de is about 2-3 (Peschard 1925). 

Neutron diffraction work on iron rich Fe—Ni alloys would provide 
a clear test of the ideas we have just described. Shull and Wilkinson’s 
experiments have been confined to alloys with less than 50%, iron and 
his results are reproduced in fig. 2 (b). The curves we would expect on 
these ideas have been roughly sketched in also. We expect jiy; to fall 
rapidly over the region of concentration where iron changes over from low 
3d behaviour to high 3d behaviour. It is difficult to know in this case 
whether the change over should be associated with the phase transition 
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or the maximum in the magnetic moment curve as these occur at only 
slightly different concentrations. The sluggish phase transformations 
would make the experiment a difficult one. 

Fe-Cu. This alloy we have included in this section because it is 
an alloy between Fe and a high 3d element. If the Cu 3d shell remains 
full, Cu will merely dilute Fe, like a Type A alloy. Unfortunately 


copper has such a low solubility in iron that it seems unlikely that any — 


experimental results will be obtained on this alloy. Gold, however, 
does act as a simple diluent (Fallot 1936). 

Fe-Zn. Similarly we expect Zn merely to dilute Fe, i.e. dji/dc=—2-2. 
Experimentally the value observed is —2-0 (Fallot 1936). 

Ni-Mn. Manganese is similar to Fe so that for small Mn concentrations 
we expect Mn to have a high-3d character. Ordered Ni, Mn is known to 
have a mean moment per atom of 1-015 (Kaya and Nakayama 1940) and 
Shull and Wilkinson (1955) has determined that either 


iy, = +3°18(+.0-25) and fiy,=-+0-30(--0-05) 
or iiwn== — 1-15(-0-25) and jiy,=1-73(+-0-05). 


The first set of results seem the more plausible ones. From these results 
Mn goes into Ni in nearly the pure configuration (3d)? with a moment 
of about 3 pp. 

The moment of the disordered alloy at room temperature initially 
rises with the addition of Mn but then passes through a maximum and 
falls rapidly to zero at about 25°%, Mn. This fall is probably partly 
due to the very rapid fall in the Curie temperature with Mn concentration 
which makes it very difficult to extrapolate to zero temperature to 
obtain the mean moment per atom at saturation. It seems natural to 
suppose that Mn enters the disordered alloy just as it enters the ordered 
alloy, as (3d)’, giving an initial value of +3-0 for dj/de consistently 
with Sadron’s result. At higher concentrations the probability of 
Mn—Mn nearest neighbours increases and the Mn—Mn coupling is probably 
antiferromagnetic. This causes the mean moment and the Curie tempera- 
ture to drop rapidly until both are zero at about a concentration of 
about 25%. In the ordered Ni,Mn alloy there are no nearest neigh- 
bour Mn—Mn pairs so the ordered alloy is ferromagnetic. An increase 
in Mn content quickly reduces the moment to zero at about 40° Mn. 

Co-Mn. Mn enters Co as it enters Ni, in the (3d)? configuration, 
but this time antiparallel to the Co moments to give a value of dji/de 
equal to about —4-0. Crangle (1957) has found experimentally a value 
of —4:5. In view of the theoretical uncertainty connected with the 
orbital contribution in the (3d)? configuration we do not think the 
difference between these results is significant. Crangle also reports 
a value of —2-0 for Mn in face-centred cubic cobalt, but this result is 
derived from a long temperature extrapolation. The reliability of such 
extrapolation in a system where Brockhouse (1957) has reported compli- 
cated antiferromagnetic coupling at low temperature is questionable. 
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An alternative explanation of the result dj/de—=—4-5 in terms of Mn 
as a low-d solute, analogous to that for Co-Cr, would give jy, 1-75, 
Moo l7—4-5c with Mn in the mean configuration (3d)!75(4ce)>, 
Neutron diffraction would readily distinguish between the suggested 
explanations. 


§ 5. CONCLUSION AND DiIscussIoN 


It is suggested in this paper that the electron distribution round 
atoms in a metal is very different from that in a free atom of the same 
element. In particular we have traced some of the consequences of 
extrapolating Weiss’s results on the number of d electrons in various 
pure metals. In conformity with his work the elements to the left 
of iron have been denoted as ‘low-d’; those to the right as ‘ high-d ’ 
configurations. The dilute binary alloy systems may be correspondingly 
grouped into four sets; low-d with low-d, and high-d with high-d, low-d 
in high-d, and high-d in low-d. The behaviour of the first two sets 
is simple, each atom retaining its own configuration in the alloy; in 
the high-d matrix some low-d elements retain their own configuration 
but Fe and Mn (which are near the dividing line between high-d and low-d) 
changeover to a high-d configuration. High-d atoms near the dividing line 
(Ni and Co) in a low-d matrix lose one or. two d electrons, the others 
(Cu and Zn) retain their usual configuration. 

To these inferences about the electronic structure of the alloys we 
have added certain assumptions concerning the magnetic moment in 
ferromagnetic alloys. These are that the magnetic moment is almost 
entirely the spin moment of electrons in 3d orbitals; that in iron the 
3d electrons are not in equilibrium with the conduction electrons, but 
that in Co and Ni there is such an equilibrium. This equilibrium in 
Ni and Co exhibits itself as a suppression of d electron spin moments 
as conduction electrons are added to the system. Conduction electron 
density changes do not influence the moment in iron. In Fe, Co and Ni 
any d electron spins carried by the solute atoms may contribute either 
ferromagnetically or antiferromagnetically to the total moment. 

With these assumptions it has proved possible to get quantitative 
agreement with experiment for all the neutron diffraction experiments 
which have been carried out on ferromagnetic alloys, as well as with 
the 22 of the 24 observed values cf dependence of saturation magnetization 
of ferromagnetic alloys. The two outstanding discrepancies are Fe-Mn 
and Co Al; both systems have maximum primary solubility ranges of 
less than 5%. 

A limited degree of success can be reported for alloys involving metals 
from the later periods ; the alloys with the Pt group of metals do not 
readily fit into the scheme ; most other available results do. 

The most controversial suggestions are that there is little occupation 
of d states until we reach Cr, and that there is a major break in electronic 
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structure after iron. This is quite different from the type of model 
in which a d-band of roughly constant band form is filled steadily from 
Ti to Cu. We make no theoretical apology for this assignment of 
structures, though it is clearly necessary to establish by calculation 
whether our opinion is in fact reasonable. A less extreme suggestion, 
but still a little surprising, is that the electronic structure is not appreciably 
affected by crystal structure. It seems to us that the nuclear charge is 
the main determining variable which, together with boundary conditions 
appropriate to a metal, fixes the electron distribution round any atom. 
The arrangement, and to some extent, the chemical nature of its neigh- 
bours are of secondary importance. The crystal structure is chosen 
to be most appropriate to the electronic structure, rather than vice versa. 
The small effect of the chemical nature of the solvent is probably a 
consequence of the rather small range of electronegativity of Fe, Co, and 
Ni. It will not be surprising if extension to other solvents requires 
somewhat greater variability in the structure of solute atoms. As an 
extreme example of this, there can be no doubt that chromous and 
chromic, ferrous and ferric ions in paramagnetic salts retain only d 
electrons. 

In adopting this model of the transition elements we lose some of 
the attractive features of the older d-band scheme. The high resistivity 
of Co and Ni compared with copper is still retained, and their high 
electronic specific heats are still given essentially the same explanation 
in terms of a high density of states in the d-band. But both resistivity 
and electronic specific heat are high in Fe, where we suppose no d-like 
high-density of states region of Fermi surface (see Mott and Stevens 1957). 
Equally the wide variation in these quantities from Ti and V to Cr is 
puzzling—as it was on the band theory. However, resistivity estimates 
for polyvalent metals are very difficult, and Jones (1957) has pointed out 
ambiguities in the interpretation of the linear specific heat term. Detailed 
work is necessary to establish whether or not any of these difficulties 
will prove unsurmountable. 

It should be pointed out that one serious objection to this model is 
concerned with the g-value of iron. Experimentally it is known that the 
spectroscopic splitting factor g is greater than two, i.e. the magneto- 
mechanical ratio is less than two. The small orbital momentum is 
aligned parallel to the spin. On our model the d-shell is less than half 
full and might be expected to lead to the opposite result. However, the 
d-shell does not now exist in isolation. The conduction electrons un- 
doubtedly carry some spin and some orbital momentum. Lack of 
detailed knowledge of the conduction band precludes detailed discussion 
of the interactions between it and the d-shell, which may invalidate the 
naive conclusion. 

Clearly, a great mass of information on transition metals and their 
alloys could profitably be re-examined in the light of these proposals. In 
addition, confirmation and extension of the Weiss ‘d electron count ’ 
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experiments would be welcome. The Fe—Co system should be of parti- 
cular interest ; according to Weiss’s scheme, the mean number of d 
electrons per atom should increase from 2 to 8 or 9 between pure Fe and 
pure Co, and in our view, this should occur precipitously between 15°/, and 
35% Co. Another type of experiment would concern the introduction 
_ of strong disorder scattering of x-rays from the widely different numbers 
of d electrons on the atoms. For example in Fe Co we suggest a difference 
of five d electrons, a very considerable difference in scattering power. 
Since the conduction electrons will also be localized to some extent to 
compensate these large differences, the effect must be looked for at 
large scattering angles where the form factors for d and s electrons are 
very different. 

Any form of test of these ideas must be welcomed, for though the 
suggested electronic structures seem at first sight quite unacceptable, 
Weiss’s results and the high degree of success of the present interpretation 
of magnetic properties are difficult to ignore. 
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Tue abundant evidence for the lack of conservation of parity in weak 
interactions leads us to ask whether parity is conserved in strong inter- 
actions. 

The spin of the “Be nucleus is 3 and the parity odd (Ajzenberg and 
Lauritsen 1955). Since the centrifugal barrier inhibits thermal neutron 
interactions with orbital angular momentum greater than zero, the parity 
of the initial system formed by the absorption of a thermal neutron must 
be odd. The symmetry and parity operators are identical for identical, 
spinless particles and so a *Be state of odd parity cannot break up into 
two alpha-particles if parity is conserved. The Q-value of the 7Be(n, «)*He 
reaction is (from the known masses) 18-98 Mev (Ajzenberg and Lauritsen 
1955) ; the disintegration of 7Be by slow neutrons into two 9-5 Mev alpha 
particles would, therefore, indicate a violation of parity conservation : 
we have sought it. 

A thin *Be target was made by evaporating down a solution of 7Be in 
HCl. The target was mounted in the dividing wall of a xenon gas 
scintillation chamber each half of which was viewed by an RCA 6342 
photo-multipler tube. A 0-02 mil nickel foil served as a light shield 
between the two halves of the chamber. The chamber contained two 
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retractable *8'Pa sources for calibration}. During the run the chamber 
was calibrated about once every two hours. Drifts were less than 
2°, per hour. 

The chamber was placed with the target above the beam hole at the 
medical facility of the BNL research reactor. The outputs of the photo- 
tubes were amplified and fed into a fast-slow coincidence circuit with the 
slow channels biased to accept all pulses above 5 Mev. The coincidence 
circuit output gated a 100 channel analyser, and also a two dimensional 
display (X-Y-Z analyser) (Grodzins 1955). Pictures on the X—Y-Z 
analyser were taken for 15-minute periods. 
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100 CHANNEL SPECTRUM 
IN COINCIDENCE WITH >5 Mev 


ye 
: \ fit 


CHANNEL NUMBER 
Spectrum observed on 100 channel analyser on one side of xenon gas scintillation 
chamber in coincidence with all pulses greater than 5 Mev on the other 
side. Bracket centres on the channel number corresponding to 9-5 Mev 
alpha particles. 

The spectrum as observed on the 100 channel analyser is shown in the 
figure ; the bracket centres on 9:5 Mev. An upper limit of about 100 
events leading toa 9:5 Mev alpha group can be set from these data. 
The large number of big pulses may be due to fission from an as yet 
undiscovered source. 


+ The authors are indebted to Dr. Elizabeth Baker and Mr. Anthony Fasano 
for making the target, and to Dr. Brian Pate for preparing the *°'Pa sources. , 
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The coordinates of all the events recorded by the X—Y—Z analyser 
were read and the results are partially summarized in the table. 

An upper limit of about 30 9-5-9-5 Mev alpha-alpha coincidence 
events can be set from these data. From the measured flux 
(5x 107 neutrons cm2/sec), 7Be target strength (4 millicuries or 10'° atoms) 
and running time (7 hours) a corresponding upper limit of 25 millibarns 
can be set for the (n, «) cross section. 

The disintegration of 7Be by slow neutrons has been previously investi- 
gated by Hanna (1955) who measured a cross section of 5 x 10* barns for 
the (n, p) reaction and who gives an upper limit of one barn for the 
(n, «) cross section. 


Total Number of Events Observed in Various Energy Intervals 


Energy (Mev) Y=5-0-6:8 6-8-8-6 8-6-L0-4 10-4-12-2 


X =6-5-11:5 85 


Following the notation of Lee and Yang (1956), the relative amplitude 
of the parity non-conserving component of the wave function is designated 
as ff. A reaction which violates a parity selection rule will have a 
relative intensity $?. So here f?~y,7/yp? where the y’s represent the 
reduced widths for the respective reactions. (We are ina region of strongly- 
overlapping resonances.) The observed width J’ is related to the reduced 
width by: I’=2kPy? where k=wave number, and P=penetration 
factor which, in both our cases, can be approximated by unity. We 
therefore find, neglecting fine points : 


2 I, 2 ky ae 
Sf Sip ies <1 mae 

Tanner (1957) puts an upper limit on sf? of ~4x10-® by failing to 
observe a resonant yield in the ground state « group at the 340 key 
J=1-+resonance in the /*F(p,«)!O reaction. Wilkinson (1958) finds 
an upper limit to $f? of ~1 x 10~7 in the *He(d, y)®Li reaction, #2<1 x 10-7 
and f?<3x10-% in the !'B(p, p’)4*B and 19F(p, «)16*O reactions 
respectively. Combining all the results to date, an upper limit of 
#?<1x 10-7 can be set for strong interactions. 

We may finally note that our result is not sensitive to whether or not 
the parity non-conserving forces conserve isotopic spin. This was not 
true of the earlier F(p, «)'®O or *He(d, y)®Li investigations into pacity 
selection rules. 
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On the Fermi Surface of Copper 
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PrpparRp (1957) has recently given, on the basis of experimental 
measurements, a description of the Fermi surface of copper. 

It is interesting that these results can be fitted to 1°/ accuracy with 
a simple analytical function. The energy £ as a function of the compo- 
nents of the electron wave-vector k must satisfy the correct boundary 
and symmetry conditions. A Fourier analysis gives terms which are 
equivalent to taking succesive orders of nearest neighbours in the 
‘tight-binding’ model, as in the interpolation scheme of Slater and 
Koster (1954). 

For a face-centred cubic crystal of lattice constant a, the first two terms 
must be of the form : 


E,=2[(—3-+cos tak, cos tak,+cos sak, cos 4ak,+cos dak, cos sak) 
+r(—3-+cos ak,+ cos ak,+ cos aks)]. 


The shape of the surface H=const. (H,,) is determined by the values of 
the parameters. We choose r=0-0995 and H/x—=3-6301, so that the radius 
vector in the (001) direction is given exactly and the experimental shape 
of the surface is fitted as closely as possible. This was verified by directly 
checking 34 experimental points distributed in such a way that, owing 
to the cubic symmetry, this amounts to studying 162 points in one 
octant. The maximum error is found to occur at the (110) direction 
and amounts to 1% of the experimental radius vector. The accuracy 
improves in going away from this point and the fit is practically perfect 
over a great deal of the solid angle, in particular in the contact around the 
(111) direction. 
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Tue efficiency of the zone refining process can obviously be increased 
by stirring the molten zone to disperse the impurity rich layer at the 
solid—liquid interface. Induction heating is sometimes preferred to 
radiant heat because it produces more convection; but no marked 
improvement has been reported. Recently Pfann and Dorsi (1957) have 
described a method of stirring the melt by passing an electric current 
through the ingot and impressing a magnetic field across the molten 
zone. The purpose of this letter is to describe the preliminary results 
obtained using a rotating magnetic field as the stirring agent, during the 
purification of aluminium. 

The metal, in the form of a rod 3in. in diameter, 44 cm long and 
99-99°%, pure, was melted in an alumina boat lined with fused alumina 
powder (which prevented the aluminium sticking to the boat). This, 
supported inside a transparent quartz tube carried by a trolley, was 
drawn through a resistance furnace. Round the furnace were mounted 
three water-cooled copper wire coils at angles of 120° to one another. 
These were fed with 50 cycle three-phase a.c. adjusted so that each gave 
approximately the same magnetic field at peak voltage. The resultant 
rotating field, of roughly 100 oersteds, was sufficient to produce visible 
stirring in mercury and to spin a large ball bearing. The axis of rotation 
was perpendicular to the rod axis so that the boundary layer was dispersed 
by the primary circulation induced in the melt. 

During operation the length of the molten zone was maintained roughly 
constant by varying the power input to the furnace, according to a 
calibration previously made on a stationary ingot. 

In order to test the effect of the stirring, the recrystallization properties 
of a rod refined with stirring were compared with those of an unstirred 
ingot. Hach rod was given eight horizontal passes in air, at a speed 
of 2:2 cm per hour; the aluminium being cleaned in aqua regia initially 
and after the 3rd, 6th and final passes. In this procedure the impurity 
rich ends were not cropped. The rods were then cut into a number of 
pieces, each of which was cold rolled 85°, with intermediate quenching 
in cold water, and etched in Lacombe’s reagent at room temperature. 
As the specimens were single crystals initially, it was not necessary to 
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polish the surfaces before etching ; but the progress of recrystallization 
could be detected by treating the material in the rolled condition. 


8 passes 


Time of recrystallization against 
distance of specimen from the end 
of the ingot 


O= Stirred ingot 


Time in hours 
WwW 


x = Unstirred ingot 


s) 10 20 30 40 L 
Pure Distance along ingot in centimetres Impure 
end end 


In the figure presented the time taken for complete recrystallization 
is plotted against the distance of the corresponding specimen from the 
end of the ingot. For the central portions of both ingots recrystallization 
was complete beween rolling and first etching (several minutes), but the 
curves for the ends show a significant increase in the efficiency of segre- 
gation of the impurities as a result of stirring. 
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REVIEWS OF BOOKS 


A History of Mechanics. By Runt Duaas, translated by J. R. Mappox. 
(London: Routledge & Kegan Paul, Ltd., 1957.) [Pp. 61] £5 Os. Od. 


THE original French edition (1950) of M. Dugas’ book has already established 
itself as a classic in the history of science, and the publication of this English 
translation is most welcome, especially as recent work on the history of mech- 
anics is too little known in this country. Before this appeared, the only treat- 
ment of the subject in English was Ernst Mach’s celebrated volume ; excellent 
as this is scientifically, and of undiminished interest for its author’s philosophical 
attitude, Mach’s history has long been out of date. M. Dugas, it is true, is not 
entirely beyond criticism on this score. The first 100 pages or so of this book 
are based on the researches of the great scholar-physicist Pierre Duhem ; but a 
great deal of further work has been done in the last forty years and Duhem’s 
conclusions have not gone unchallenged. It is a pity also that the illuminating 
studies of M. Koyré (1939) were not more fully used. Nevertheless, it is 
extremely gratifying to have a coherent survey of statistics and dynamics 
down to the time of Galileo that is representative of modern knowledge : this 
has not been available in English before. This early story occupies 128 pages ; 
330 are devoted to classical mechanics; and the final 178 pages to ‘ The 
Principles of Modern Dynamics’. Thus, as seems appropriate, the greater 
emphasis of this History is on the classical theories, but the twentieth century 
is also treated in considerable detail. Though some continuity in tracing the 
development of concepts or techniques is lost in a biographical arrangement, 
M. Dugas’ exposition is always clear, and well supported by quotation from the 
originals. The translation is in general adequate; the rendering of some 
proper names (Cues for Cusa ; Erasme for Erasmus ; Huyghens, an erroneous 
spelling), and the capitalization is not in accord with normal English usage. 
Apart from a rather unnecessary frequency of misprints the Swiss printers have 
produced an agreeable volume. It is one which everyone with an interest in 
the history of science will wish to read, and consult ; and which may well prove 
to have more than occasional value in the teaching of mechanics. 
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ERRATA 


The Surface Energy of Colloidal Metals in Ionic Lattices, by A. B. Scorr, 
1954, Phil. Mag. (7), 45, 610. 


In § 5, the potential energy of two square 36-ion planes, separated by 
ro, Should be divided by 2 before comparison with the energy of an ion- 
plane with its image plane, since the separation of the two ion-planes 
creates twice the crystal surface that is formed when an ion-plane is 
separated from the metal. The image effect in table 1 should be every- 
where doubled and the energy per atom is : 


The observed heat of dissociation then corresponds to a particle of 
about 42 atoms of K and 20 atoms of Na (page 617). 

Page 616, last line, for 294 read 296. 

Page 617, first line, for 84 read 80. 

Page 618, 18th line, for 0-001 ev read 0-06 ev. 

Figure 1, ordinates are in units of 8 e?/ro. 


I am indebted to Professor F. G- Fumi for calling my attention to these 
errors. 


[The Editors do not hold themselves responsible for the views 
expressed by thew correspondents. | 


